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ABSTRACT 

Context. We report on the results of the first XMM-Newton systematic "excess variance" study of all the radio quiet. X-ray un- 
obscured AGN. The entire sample consist of 161 sources observed by XMM-Newton for more than 10 ks in pointed observations, 
which is the largest sample used so far to study AGN X-ray variability on time scales less than a day. 

Aims. Recently it has been suggested that the same engine might be at work in the core of every Black Hole (BH) accreting object. 
In this hypothesis, the same variability should be observed in all AGN, once rescaled by the Mbh (Mbh) and accretion rate (in). 
Methods. We systematically compute the excess variance for all AGN, on different time-scales (10, 20, 40 and 80 ks) and in different 
energy bands (0.3-0.7, 0.7-2 and 2-10 keV). 

Results. We observe a highly significant and tight (~ 0.7 dex) correlation between cr^^^^ and Mbh - The subsample of reverberation 
mapped AGN shows an even smaller scatter (only a factor of 2-3) comparable to the one induced by the Mbh uncertainties. This 
implies that X-ray variability can be used as an accurate tool to measure Mbh and this method is more accurate than the ones based 
on single epoch optical spectra. This allows us to measure Mbh for 65 AGN and estimate lower limits for the remaining 96 AGN. 
On the other hand, the cr^,^j vs. accretion rate dependence is weaker than expected based on the PSD break frequency scaling. This 
strongly suggests that both the PSD high frequency break and the normalisation depend on accretion rate in such a way that they 
almost completely counterbalance each other (PSDa„,p oc m"" **). A highly significant correlation between cr^^^^ and 2-10 keV spectral 
index is observed. The highly significant correlations between cr^^^^ and both the Lboi and the FWHMh/j are consistent with being just 
by-products of the (T,-,,,^ vs. Mbh relation. The soft and medium cr^^^^ is very well correlated with the hard cr^„^^, with no deviations 
from a linear one to one correlation. This suggests that the additional soft components (i.e. soft excess, warm absorber) add a minor 
contribution to the total variability. Once the variability is rescaled for Mbh and m, no significant difference between narrow-line and 
broad-line Seyfert I is observed. 

Conclusions. The results are in agreement with a picture where, to first approximation, all local AGN have the same variability 
properties once rescaled for Mbh and in. 

Key words. Galaxies; active - Galaxies: Seyfert - quasars: general - X-rays: general 

1 . Introduction variations, i.e. tiieir PSD were well fitted by power-laws with 

slopes steeper than 1 . They also found that the PSD amplitude at 

Rapid variabdity is one of the major observational properties ^ ^^^^ frequency was anti-coiTelated with the source luminosity, 

of accretion onto Black Holes (BH). Since the early observa- confirming previous results 
tions, it was realised that the fast variability observed in AGN 

H I implied that the primary source had to be both very compact The availability of well sampled, very long (up to ~ 10-15 

- - • and able to emit with a large luminosity, thus requiring an ex- years in some cases) RXTE light curves, and their use together 

tremely efficient engine, providing the first compelling argument with high signal to noise, continuous XMM-Newton light curves 

for the presence of BH in the heart of AGN (e.g. Lynden-Bell has revolutionised our view of AGN X-ray variabiUty in the last 

1969; Rees 1984). Eai-ly studies of the AGN X-ray variabil- 15 years. For a dozen objects, the combined use of these Hght 

ity properties suggested that more luminous sources are "less" curves allowed the accurate determination of their PSD over a 

variable (Barr & Mushotzky 1986). Later on, the EXOSAT long broad range of time scales. The past suggestions that the steep 

looks allowed the first study of the observed variations with the PSD at high frequencies had to flatten below a certain time scale 

use of Fourier analysis methods, like the power specti-al density (see e.g. Papadakis & McHardy, 1995) were soon confirmed: 

function (PSD), for the brightest AGN. Lawrence & Papadakis all the best quality PSD are best fitted by a steep power law of 

(1993) and Green, McHardy & Lehto (1993) found that the ma- slope ~ -2 which "breaks" to a flatter (~ -1) slope at frequen- 

jority of the EXOSAT AGN light curves exhibited "red noise" cies below the so-called "break frequency" Vb (e.g. Edelson & 

Nandra, 1999; Utdey et al. 2002; Markowitz et al. 2004; Kelly 

Send offprint requests to: Gabriele Ponti et al. 2010). In some AGN the -1 slope part could be measured 

e-mail: ponti@iasfbo . inaf . it for more than 3 decades in frequency, with no lower frequency 
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breaks being detected, thus indicating a behaviour similar to 
Galactic X-ray black hole binary candidates (BHB) in their so- 
called "soft state" (e.g. McHardy et al, 2004). In at least one 
case, namely Ark 564, two frequency breaks have been detected 
(Papadakis et al 2002). In fact, the PSD of this object may be 
better fitted by "Lorentzian" functions (as opposed to "broken" 
power-laws; McHardy et al. 2007), similar to the functions that 
are widely used in the parametrisation of the BHB PSD. 

The detailed PSD analysis of the high quality RXTE and 
XMM fight curves revealed a remarkable similarity between tfie 
AGN and BHB power- spectra, although the time scales sampled 
in the two classes of objects were vastly different. This result 
implied a deep link between all BH accreting objects, reinforc- 
ing initial suggestions that both systems host the same engine, 
and that the same emission and variability mechanisms operate 
in them. The difference in time-scales could be explained by the 
difference in the mass of the compact object in the center of 
these objects. The results from these first studies made it clear 
that the observed AGN X-ray variations are most probably de- 
termined by BH mass (Mbh) rather than luminosity (Hayashida 
et al. 1998; Czemy et al. 2001; Uttley et al. 2002), and even sug- 
gested that X-ray variabifity might be used as a tool to estimate 
Mbh in these objects. 

McHardy et al. (2006) were the first to demonstrate clearly 
that AGN are scaled-up versions of BHB. Using PSD results for 
almost a dozen objects, spanning a range of ~ 8 orders of mag- 
nitude in Mbh and ~ 3 orders of magnitude in accretion rate 
(hereinafter we will use the term accretion rate as a sinonimus 
of Eddington ratio), they demonstrated that the PSD "break time 
scales" increase proportionally with Mbh. and decrease with in- 
creasing accretion rate. Koerding et al. (2007) showed that the 
three physical parameters, namely Mbh, accretion rate and break 
frequency, are intimately related, both in AGN and BHB in the 
soft state, determining a "variabifity plane" of accreting BHs. 
Interestingly, if a constant offset is introduced, this plane can 
also be extended to the "hard state" objects as well. Casella et al. 
(2008) even used this "variability plane" to estimate the Mbh of 
Ultra Luminous X-ray sources (ULXs). 

Although the PSD analysis is a powerful tool to characterise 
the variability properties of AGN and BHB, it nevertheless re- 
quires long, uninterrupted observations from especially tailored 
morutoring campaigns to fully exploit its potential. Such ob- 
servations are available at the moment for no more than two 
dozen AGN. However, the last few years, shorter, high signal- 
to-noise X-ray observations for tens of objects have populated 
the archive of X-ray observatories such as Chandra, and, per- 
haps even more importantly XMM-Newton (due to its high sen- 
sitivity, and broad energy band pass). A conveiuent analysis tool 
for such short data sets is the so-called "normalized excess vari- 
ance", CTjjjjj (Nandra et al. 1997). Although it does not offer the 
same wealth of information fike the PSD analysis, it can cer- 
tainly be used to confirm the PSD results using larger data sam- 
ples. In fact, the availability of larger samples can also allow the 
discovery of new correlations between the X-ray variability am- 
plitude and other AGN physical parameters. 

The first AGN excess variance surveys were performed using 
ASCA fight curves. The results confirmed the variabifity vs. lu- 
minosity anti-correlation (Nandra et al. 1997; George et al. 2000) 
but also indicated, for the first time, that the variability amplitude 
(i.e. the excess variance) correlates with the X-ray spectral in- 
dex, and anti-correlates with the FWHM of the HjS line (Turner 
et al. 1999). Similar studies also indicated, for the first time, that 
the so called "Narrow Line Seyfert 1" galaxies were systemati- 
cally "more variable" than the classical AGN of equal luminos- 



ity (Leighly 1999). Later on, the variability vs. luminosity anti- 
correlation was also confirmed on long time scales (Markowitz 
& Edelson 2001), and it was soon suggested that the excess vari- 
ance vs. luminosity relation might be just a by-product of a more 
"fundamental" relation, that of the "variability vs. Mbh relation" 
(Lu & Yu, 2001; Bian & Zhao, 2003; Papadakis, 2004). These 
suggestions were put forward even before the McHardy et al 
(2006) and Koerding et al (2007) Mbh (and accretion rate) scal- 
ing relations were published. These relations, together with the 
standard accretion theory of a-disks (Shakura & Sunyaev 1973), 
which predicts that all the disc characteristic time scales should 
depend linearly on Mbh (e.g. Treves, Maraschi & Abramowicz 
1988), reinforce the hypothesis of the "variability vs. Mbh" re- 
lation being the fundamental relation which determines all the 
other observed relations, although alternative suggestions are 
still being considered (Liu & Zhang 2008). 

The largest collection so far of excess variance measure- 
ments in local AGN has been presented by O'Neill et al. (2005). 
These authors used aU the available 40 ks long A5CA light curves 
to calculate the excess variance of 46 AGN. They found a strong 
anti-correlation of variability with Mbh and with luminosity. 
However, the latter correlation disappeared once the cr^^^ vs. 
Mbh correlation was taken into account. O'Neill et al. (2005) 
found a weaker correlation (~ 96 %) between cr^^ and 2-10 
keV spectral index, F, than observed previously (Turner et al. 
1999). Miniutti et al. (2009) confirmed the variability vs. Mbh 
trend extending the variability estimation to smaller Mbh AGN. 
Finally, Zhou et al. (2010) tried to cafibrate accurately, and esti- 
mate the intrinsic scatter of the cr^^ vs. Mbh relation, using high 
quality XMM-Newton fight curves of AGN with Mbh measured 
through the "reverberation mapping" technique. They found that 
the intrinsic scatter of this relation is even smaller than the one 
impfied by the Mbh uncertainties. Moreover, no dependence of 
the variability on either accretion rate or spectral index was ob- 
served. 

In this paper we present the results from the excess variance 
measurements of a sample of 161 AGN. This is about three times 
larger than the ASCA sample of O'NeiU et al. (2005) and, at 
present, the largest sample in which the short time variabifity 
(less than a day) has been systematically investigated. In addi- 
tion, we estimated the excess variance in various energy bands, 
and on different time scales. The first aim of our work is to inves- 
tigate whether the results obtained from the detailed PSD anal- 
ysis of good quality light curves of a few AGN are applicable 
to the "majority" of the X-ray studied AGN as well. Note that 
the McHardy et al (2006) and Koerding et al (2007) scaling re- 
lations were based on the accurate PSD results for a relatively 
smaU number of objects (for example there were just 10 AGN in 
the McHardy et al. 2006 sample). The excess variance, being the 
integral of the PSD over the frequency window sampled in the 
fight curve that is used to compute it, can be a powerful tool that 
can be used to investigate the appficability of these relations to 
a much larger sample of AGN. This can be done with the excess 
variance vs. Mbh and excess variance vs. accretion rate relations, 
which we study in detail in this work, putting particular empha- 
sis on the comparison between the PSD model predictions and 
the observed relations. 

The second major aim of our work is to investigate how ac- 
curately we can "weigh" the central Mbh in AGN with cr^^,, 
and we provide accurate recipes to measure Mbh from X-ray 
variabifity using light curves of various lengths. Moreover, we 
present the results from the study of the variability amplitude 
with the X-ray spectral index and FWHM H/3. Our results indi- 
cate that the same variability mechanism operates in all AGN, 
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and that significant differences in the variability amplitude of 
NLSl and "typical" broad line Seyferts (BLSl) can be fully un- 
derstood once the observed variability is properly "normalized" 
to the Mbh and accretion rate of these objects. 

2. The catalogue: CAIXAvar 

XMM-Newton provides high statistics, low background, uninter- 
rupted light-curves, allowing the calculation of the excess vari- 
ance in large samples of AGN. We study here a new catalogue, 
CAIXAvar, that is a sub-catalogue of the CAIXA sample pre- 
sented by Bianchi et al. 2009a,b. CAIXA consists of all the 
radio-quiet X-ray unobscured (Nh < 2 x 10^^ cm"^) AGN ob- 
served by XMM-Newton in targeted observations. We selected 
a sub-catalogue of CAIXA, creating CAIXAvar, which includes 
datasets of sources with the following characteristics: 1) cleaned 
exposure times larger than 10 ks; and 2) at least 20 counts in the 
(rest-frame) 2-10 keV band for each time bin of 250 s. Compared 
to CAIXA (Bianchi et al. 2009a,b), we expand the sample in- 
cluding all observations whose data are public as of June 2010 
and consider here also multiple observations of the same object. 
There are 161 sources in CAIXAvar that fulfil these criteria (with 
260 observations), 125 and 158 of which have a measurement of 
the Mbh and FWHM of Hj}, respectively. All sources have at 
least one observation with exposure longer than 10 ks, and there 
are 89, 56 and 32 sources with at least one observation with ex- 
posure longer than 20, 40 and 80 ks, respectively. 

We estimate the bolometric luminosities (Lboi) through four 
different recipes. First, we collected the values reported by Woo 
& Urry 2002, whose measure of the LboI is derived by integrat- 
ing the source spectral energy distribution. The other three meth- 
ods use the mean X-ray luminosity of the longest XMM-Newton 
observation. We used the constant X-ray bolometric correction 
by Elvis et al. (1994), the luminosity-dependent one provided 
by Marconi et al. (2004) and, finally, the BH-mass dependent 
bolometric correction by Vasudevan et al. 2007. We checked and 
confirmed that the results presented here are not affected by the 
bolometric correction we use. Thus, in this work we report the 
results when we used the Marconi et al. (2004) bolometric cor- 
rection to estimate LboI for all objects in the sample. This cor- 
rection is applicable to all the sources in CAIXAvar, allowing us 
to use the same recipe uniformly throughout the sample. We also 
use these values to compute the ratio LBoi/LEdd (for the sources 
with available Mbh estimates) which we will refer to as the "ac- 
cretion rate" (m) of the source. 

Mbh from reverberation mapping are preferentially used, 
then stellar velocity dispersion, measurements using the relation 
between broad-line region radius and optical luminosity (primar- 
ily based on Hj3, but also on Mg II whenever the former was not 
available) and finally the relations for the narrow line region. 
The measurements of the stellar velocity dispersion are col- 
lected from the HyperLeda database {http://leda.univ-lyonl.fr/) 
and the Mbh estimated using the relation: LogiMsH/Mo) - 
8.12 + 4.24 Log(o-J200 km i"'); see Gultekin et al. (2009). 
However, as Greene & Ho (2006) and Greene et al. (2008) sug- 
gest, the Mbh vs. stellar velocity dispersion might deviate at low 
Mbh from this relation, thus for Log(MBH)<6 the values are di- 
rectly taken from Greene & Ho (2006). 

Table [T] summarises the multi-wavelength data collected for 
the sources in CAIXAvar. The left panel of Figure [T] shows the 
redshift distribution of all the 161 AGN in the CAIXAvai" sam- 
ple. The vast majority of these AGN are local with redshift lower 
than 0.2. However, CAIXAvar contains also 20 AGN with red- 
shift higher than 1 and in particular 4 AGN with z as high as 4. 




-1.5 -1 -0.5 0.5 -3 -2.5 -2 -1.5 -1 -0. 

Log(redshilt) Log(redshift) 



Fig. 1. (Left panel) Redshift distribution of the all the 161 
sources in CAIXAvar. Although the vast majority of AGN be- 
long to the local Universe z< 0.2, 20 sources in CAIXAvar have 
z> 1, with 4 objects with redshifts as high as 4. (Right panel) 
Redshift distributions of all the sources with at least 1 variable 
segment (see bold values in Tab. 2) used hereinafter to study the 
correlations. Only 3 objects, with at least a variable segment, 
have a redshift higher than 0.2 and only one have z ^ 0.9. 



The right panel of Fig.[T|shows the distribution of the sources of 
CAIXAvar with at least one variable segment. The sources with 
at least one variable segment (which we use to study the vari- 
ability and correlations of the CAIXAvar sample with different 
physical parameters) is primarily made of local AGN. In fact, 
there are only 3 objects with redshift higher than 0.2 and only 
one with z > 0.9 with at least one variable segment. 

2. 1. Reverberation sample (Rev) 

In order to keep the scatter on the relations with the Mbh as small 
as possible we also selected a side sample ("Rev", hereafter) of 
AGN with Mbh measured through reverberation only. Peterson 
et al. (2004) measured the Mbh for 37 AGN, 8 of which have 
either: no; or still proprietary; or too short; XMM-Newton obser- 
vations. The remaining 29 sources are analysed regardless of the 
radio loudness and absorption (after individually checking that 
these do not play a major contribution to the measured 2-10 keV 
variability). All the reverberation Mbh are from Peterson et al. 
(2004) apart from PG1351H-442, whose Mbh estimate is taken 
from Kaspi et al. (2000). The Mbh estimate for the 6 objects 
(one of which, Mrk290, is not in the Peterson et al. 2004 sam- 
ple) re-analysed by Denney et al. (2010) is taken from that paper 
Finally, the Mbh estimate for Mrk766 (which is not included in 
Peterson et al. 2004) is taken from Bentz et al. (2009). 

Due to the poor quality of the data, Peterson et al. (2004) 
mark as less accurate the measurements of the Mbh of IC4329A 
and PG121 lH-143. Following Markowitz et al. (2009) we use for 
IC4329A a value of Mbh = 2.17 x 10^ Mq, obtained through 
the stellar velocity dispersion (cr, = 218 - 231 km s"'; Oliva et 
al. 1999) that is in agreement with all the other methods used to 
estimate the Mbh- PG1211h-143, instead, has no stellar velocity 
dispersion measurement. We note that this value is almost one 
order of magnitude higher than the reverberation mapping esti- 
mate given by Kaspi et al. (2000; Log(MB//) = 7.37) and more 
than one order of magnitude higher compared to the Mbh esti- 
mates from single epoch spectra (Hyj; Log(MBH) = 6.72 + 0.43; 
CIV 6.95 + 0.40 and from BLR estimate 6.84 + 0.29; Kelly & 
Bechtold 2007). In the end, we decided to use the Peterson et 
al. (2004) measurement. We also added NGC4395 in the Rev 
sample because an HST campaign allowed the reverberation 
mapping measurement of Mbh = 3.6 x IO^Mq for this object 
(Peterson et al. 2005). For conformity we use this value even if 



3 



G. Pond et al.: CAIXA: a catalogue of AGN in the XMM-Newton archive III. Excess Variance Analysis 



Sample 


Name 


Z Log(MBH) 


Log(LBol,Woo) 


Log(LBol,Elv) 


Log(LBol,Mar) 


Log(LBol,Vas) 


FWHM 


r 






(Mo) 


(erg s"') 


(erg s"') 


(erg s"') 


(erg s"') 


(km s-') 




R+C 


MRK335 


0.0258 7.15±0.11 


44.69 


44.3561 


44.0919 


44.5322 


1620 


1.53 



Table 1. List of all Multi-wavelength data for each source. The first column indicates to which samples the source belongs (R for 
the Rev sample and C for the CAIXAvar sample). The entire table is published in the on line version of the paper, only. 



we note that it might be an overestimate of the true mass. In fact 
the low bulge velocity dispersion (cr, < 30 km s"'; Filippenko 
& Ho 2003) deviates from the extrapolation of the M-cr rela- 
tion when this estimate of the Mbh is adopted (but see Greene & 
Ho 2006). The final Rev sample consists of 32 AGN (11 more 
AGN than in the Zhou et al. 2010 sample), 6 of which are ei- 
ther absorbed or radio loud objects, hence they are not present in 
CAIXAvar 



3. Observations and data reduction 

The analysis products have been obtained starting from the ODF 
files and after reprocessing with SASv6.9. The source and back- 
ground regions have been selected with the same procedure as 
Bianchi et al. 2009a,b apart from the observations in Small 
Window mode for which the background has been selected from 
a source-free region on the same chip as the source (source and 
background have the same extraction radius). The screening for 
flaring of particle background was performed in two steps. The 
first cut has been performed in the same way as Bianchi et al. 
(2009a,b), via an iterative process that leads to a maximisation 
of the signal-to-noise (see Piconcelli et al. 2005). We also filtered 
out every timebin during which the 10-15 keV light curve of the 
entire field of view (but the source region) had more than 30 and 
1250 counts during the 250 s bins for the Small Window and 
Full Frame mode, respectively. This prevents small fluctuations 
of the particle background from making an appreciable contribu- 
tion to the excess variance measurement. Objects requiring long 
cuts have been individually screened, and an ad hoc background 
selection has been performed on a case-by-case basis. 

Light curves have been constructed in the energy bands: 
0.3-0.7, 0.7-2, 2-10 and 0.3-10 keV, by selecting events with 
FLAG==0, PATTERN<=4 and XMMEA.EP . The energy band 
limits refer to rest-frame energies (thus we were unable to com- 
pute the low energy excess variances for the 15 objects with a 
redshift larger than about 2). Light curves in all bands are cor- 
rected with the SAS task EPICLCCORR. Every light curve was 
calculated with 250 s bins, and divided into segments of 10, 20, 
40 or 80 ks. During the estimation of the excess variance (see 
below), we rejected the time bins with fractional exposure lower 
than 0.35, regardless of the observation mode (in Small Window 
mode the EPIC pn camera is active for only about ~ 70 % of the 
time, thus we reject only time bins with background flare lasting 
more than half of the active time). 

3. 1. Normalised excess variance computation 

Following Nandra et al. (1997), Turner et al. (1999), Vaughan 
et al. (2003) and Ponti et al. (2004), we compute, for each light 
curve segment, the normalised excess variance with the formula: 

1 " 



where is the number of good time bins in the segment, /i 
is the unweighted arithmetic mean of the counting rates within 
the segment, X, and cr, are the counting rates and uncertain- 
ties, respectively, in each bin (hereinafter, we refer to cr^^^^ sim- 
ply as "excess variance" of the light curve). We rejected all the 
segments with less than 50 % of good time bins. When more 
than one valid segment is available, the excess variance has been 
determined computing the unweighted mean of all the individ- 
ual estimates. This has the potential to reduce the large uncer- 
tainty owing to the stochastic nature of the excess variance. The 
method of the excess variance uncertainty estimation and the 
impact of different source redshifts is described in detail in the 
Appendix A and B. 

4. The excess variances 

Table 12] shows the excess variance estimates and the source pa- 
rameters for the CAIXAvar and Rev samples, respectively, in the 
2-10 keV band. 

In several cases, due to the fact that we were quite conser- 
vative in assigning the 90% confidence limits on our estimates, 
and especially for shorter intervals (10-20 ks) and less variable 
AGN (which are those objects with the largest Mbh), the lower 
limit on cr^^^^ implies an intrinsic excess variance less than zero. 
In this case, we consider our measurement as a "non-detection", 
and we simply list in the respective tables the 90% upper limit 
of our estimate. 

4. 1. Metliod of analysis 

In the following sections we study the correlations between the 
excess variance in various energy bands, and also between cr^^^ 
and other source parameters. We used various methods (such as 
a linear bisector and/or a Y/X linear model fit to the data in the 
log-log space, methods based on the points scatter, etc) to mea- 
sure the correlations between two parameters. We observed that, 
despite some slight differences, all the methods yielded simi- 
lar results. We thus decided to apply a censored fit of a linear 
model with the bisector method to all the data sets we consid- 
ered (see Appendix A of Bianchi et al. 2009), mainly for con- 
sistency with previous results (e.g. Bianchi et al. 2009a,b). All 
model fits were performed in the log-log space, and in all cases, 
the excess variance is considered to be the "dependent" vari- 
able, and the other parameters the "independent" variable. We 
performed MonteCarlo simulations to compute the uncertainties 
on the best-fit parameters and the correlation significance as in 
Bianchi et al. (2009). 

Finally, following Nandra et al. (1997) and O'Neill et al. 
(2005), in our model fits we considered only the excess vari- 
ances of the objects with at least one segment with significant 
variations (we also plot only these objects in all figures below). 
To identify these "variable" segments we performed ax^ test to 
each one of them (for each source and any interval length, 10, 20, 
40 and 80 ks), accepting as evidence for "significant" variabil- 
ity a null hypothesis probability of less than 1%. We measured 
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Sample 


Name 


Log(MBH) 

(Mo) 


Log(LBoiMar) 

(erg s"') 


^ mis,80,2-10keV 


" rms,40.2-IOkeV 


^ rms,20,2-10kcV 


rms,10,2-10keV 


R+C 


MRK335 


7.15+0.12 


44.09 


2.6!/,-^ X 10-^ 


2.3+/|-^ X 10-^ 


1.5 ± 0.2 X 10-- 


1.02 ± 0.16 X 10-- 



Table 2. List of all the cr^^ 2-iokeV computed, in the 2-10 keV band, with 10, 20, 40 and 80 ks intervals. This is an update of the 
multi-wavelength table published in Bianchi et al. (2009). The entire table is published in the on line version of the paper, only. 



at least a "variable" interval in 25, 42, 54 and 47 sources in the 
80, 40, 20 and 10 ks intervals, respectively. In total, there are 65 
different sources (i.e. 40 % of the 161 sources in the CAIXAvar 
sample) with at least one variable segment in at least one time in- 
terval. Table 2 lists in boldface the crj^^^ measurements for these 
sources, which we use in the studies we present here. For the re- 
maining 96 sources we do not detect significant variability dur- 
ing any segment and the corresponding upper limits are reported 
in Table 2. Not surprisingly, it turned out that most objects with 
upper limits on cr^^^ (in the case of "non-detections") did not 
have a single "variable" segment. This fact implies that we did 
not consider most of the objects with upper limits on their cr^^^ 
in the model-fits. This is not a serious problem as in most cases 
the upper limits are so large that their presence in the correlation 
plot do not add any physically constraining power. Nevertheless, 
we always checked "a posteriori" that none of these discarded 
upper limits is observed below or close to the best fit relations 
(thus carrying potentially important information). 

5. Comparison of o-^j, in various energy bands 

The three energy bands (0.3-0.7, 0.7-2 and 2-10 keV) we con- 
sider in this work are dominated by different spectral com- 
ponents (i.e. warm absorber and the so-called soft excess are 
more pronounced in the softer energy bands). Their presence 
can introduce additional spectral variability in these bands. 
Consequently, a non-linear relation between the measured crf^^ 
in these bands may not be surprising. Even if the best-fit relation 
is linear when we consider the cr^^^^ in one band vs. cr^j-,,, in an- 
other, a departure from a slope of ~ 1, or a best-fit normalisation 
different than unity, could be expected. 

The upper left and right panels of Figure |2] show a plot of 
'^rms,o.3-o.7keV f^Ls,2-iokeV (^r the CAIXAvar sample), when 
computed within time intervals of 80 ks and 40 ks, respectively. 
Similarly, the lower panels show the same plots for the 0.7-2 vs. 
2-10 keV bands. Clearly the excess variance in the soft bands are 
well correlated with cr^ t .at ^;■ Our best-fit results are listed in 

rm.s,2-10keV 

Tab.|5] The best-fit models are also plotted in Fig.|2](solid lines). 
The dotted lines represent the combined 1-cr error on the best-fit 
slope and normalisation values, while the red dashed line shows 
the one to one relationship. 

The best-fit slope and normalisation estimates are consistent, 
within the errors, with 1 and 0, respectively, in the case of the 
excess variance estimates which are calculated using the 40, 20 
and 10 ks light curves. This result implies that, on these "short" 
time scales, the main driver for the X-ray variability in these ob- 
jects is the continuum normalisation variations, while variations 
of the other spectral components must be of a much smaller am- 
plitude. In the case of the 80 ks based cr^^^ estimates, the best-fit 
results suggest a slope steeper than 1 and a normalisation larger 
than zero (in the log-log space). A steeper slope implies that the 
more variable AGN tend to sit above the 1 to 1 relation, while 
less variable objects below. Consequently, the PSD shape and/or 
amplitude may not scale from one energy band to the other in 
the same way for all objects and, as a consequence, we observe 



a not null normalisation once the relation is extrapolated to 0. 
However, the significance of this result is rather low (at the ~ 2cr 
level) . 

The very good correlation between the excess variance mea- 
surements in the three bands could allow us to integrate the sig- 
nal over the largest possible energy band (i.e. 0.3-10 keV) in 
order to gain better precision on the variability measurements. 
However we present below the results from the correlation of the 
2-10 keV band excess variance measurements with other source 
parameters, for comparison with previous work, and in order to 
allow a future extension of this work using longer timescales 
data from the RXTE PCA observations. 



6. Correlations between thie excess variance and 
othier source parameters 

To investigate potential physical mechanisms which drive the X- 
ray variability in local AGN on short time scales (i.e. less than 
a day) we search for correlations between cr^^ (i.e. the source 
variability amplitude) of the variable sources in CAIXAvar and 
various physical parameters of the sources such as: Mbh, accre- 
tion rate, luminosity, width of the H^ line, and X-ray spectral 
slope. We use the excess variance measurements from the 10, 
20, 40 and 80 ks segments when we study the variability vs. 
Mbh or m relation. In all other cases, for brevity reasons, we use 
the cr^jj^j, measurements from the 20 and/or 40 ks intervals only. 

CAIXAvar is the largest sample used so far to study the AGN 
X-ray variability on short time scales (less than a day). Large 
samples of AGN have been used to study the variability prop- 
erties of higher redshift AGN, and on longer time scales, albeit 
with light curves which are heavily undersampled. For exam- 
ple, Almaini et al. (2000) and Manners et al. (2002) have used 
ROSAT data to study the variability properties of 86 and 156 
AGN, respectively. Both Paolillo et al. (2004) and Papadakis et 
al. (2008), thanks to the study of variability of Chandra and 
XMM-Newton deep field observations, trace a higher redshift 
AGN population. Paolillo et al. (2004) study a sample of 430 
Chandra sources from the deep field south (74 of which result 
to be variable) and Papadakis et al. (2008) 66 AGN from the 
XMM-Newton Lockman Hole observations. Both these works 
study AGN variability, but only on time scales longer than ~ 1 
day. Finally, Vagnetti et al. (2011) used a different technique, 
the structure function, to study the variability of a large sam- 
ple (412) of AGN from the XMM-Newton Serendipitous Source 
Catalogue. 

6.1. Thecrl^^ vs. Mbh relation 

The presence of a correlation between X-ray variability and Mbh 
is ah-eady well established (Lu & Yu 2001; Bian & Zhao 2003; 
Papadakis 2004; O'Neill et al. 2005; Nikolajuket al. 2006; 2009; 
Zhou et al. 2007; Miniutti et al. 2009; Zhou et al. 2010). First, 
we study this correlation using the results of the Rev sample, for 
which every object has a well determined Mbh- This allows us 
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o'LCS-IO) 

log((7^„3(07-2) ) = (0.42±0.19) + ( 1.19±0.10) log((T^„3(2-10)) [OLSB] 
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Fig. 2. Upper panels Soft (0.3-0.7 keV) vs. hard (2-10 keV) cr^^^, computed within the 80 ks (left), and 40 ks (right) intervals. 
Lower panels. Similar plots for cr^^^ in the 0.7-2 vs. cr^^^ in the 2-10 keV band. The best fit curves are plotted with solid lines and 
the combined 1-cr error on the slope and normalisation with dotted lines. The red dashed lines represents the one to one relations 
expected in the case of achromatic variations. On top of each panel (in this and all similar subsequent figures) we report the best fit 
models (see also Tab.|5]l. 



to keep under control the scatter introduced by the mass uncer- 
tainties, that are estimated to be of the order of a factor of 3 or 
smaller, for this method (Peterson et al. 2004). Then, we extend 
our study to the full CAIXAvar sample, using the sources in this 
sample with a Mbh estimate. 



6.1 .1 . Results on the Rev sample 

The four panels of Figure |3] show the cr^nis 40 20 lo ^bh 
plot for the Rev sample, respectively (cr^^, o"?ms,40' ^rms,20 ^^'^ 
•^rms 10 denote the excess variance measurements from the 80, 
40, 20 and 10 ks segments, respectively). A highly significant 
anti-correlation is observed in all cases (see Tab.|3]l. Solid lines in 
these panels indicate the best-fit relations, and the best-fit results 
are listed in Tab. [3] We note that, as expected, the normalisation 
of the relation increases for longer intervals. As the light curve 
segment duration increases, the expected variance is expected to 
increase as well, as in effect it is like integrating the PSD over 
a larger frequency window, hence we do expect to measure a 
larger variability amplitude. 



To quantify the scatter of the data around the best- 
fit line we compute the quantity (Xscatter as follows: 

cTscatter = ^ZliiLogicT^^^i) - Log[f(M bhjW /N; where 
cr^^_,j and /{Mbhj) are the measured and best-fit excess vari- 
ances for a source with Mbh of Mbh, respectively, and is the 
total number of sources, excluding sources with upper limits on 
their cr^^^. We also exclude PG1211H-143; this source is the only 
clear outlier in the middle and right panels of Fig. 3, but also has 
the largest uncertainty associated with its Mbh measurement 
(see ^T[ . 



The o-scattei-,80,40 and crseatter,20 values are 0.44, 0.49 and 0.47, 
respectively, and they imply a scatter of a factor of 3 in linear 
space. Interestingly this is comparable with the uncertainty asso- 
ciated with the reverberation mapping Mbh measurements. This 
result implies that the variability Mbh estimates, using the best- 
fit lines with the parameters listed in Tab.[3j are at least as accu- 
rate as the reverberation mapping ones. 
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log(<TL) = (-19''±0.12) + (-1.15±0.12) log(MB„,) [OLSB] log(a^„J = (-2.00±0.13) + (-1.32±0.14) log(Mg„^) [OLSB] 
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l°g(gL) = (-213±0.14) + (-1.24±0.12)log(M|3H^) tOLS B] log(g'^,) = (-2.20±0.14) + (-1.21 ±0.10) logiMeH,) lOLS B] 




10" 10' 10° 10° 10° 10' 10° 10' 



Fig. 3. Upper left, upper right, lower left and lower right panels show the o-^^^^q, o"rms40' "^^520 ''"rms 10 for the Rev 

sample. The best fit relationships (see Tab.O are plotted with solid lines and the combined 1-cr error on the slope and normalisation 
with dotted lines. Mbhj indicates Mbh in units of 10^ M©. 



Relation 


80 ks 

Norm Slope prob 


40 ks 

Norm Slope prob 


Rev 

cr;^, vs. Mbh 


-1.94 + 0.12 -1.15 ±0.12 99.998 


-2.00 ±0.13 -1.32 ±0.14 99.998 


Relation 


20 ks 

Norm Slope prob 


10 ks 

Norm Slope prob 


Rev 

f^L vs. Mbh 


-2.13 ±0.14 -1.24 ±0.12 99.79 


-2.20 ±0.14 -1.21 ±0.10 99.6 



Table 3. List of all best fit relations of the Rev sample as well as their probabilities. 



6.1 .2. The full sample: CAIXAvar 

The four panels of Fig. g] show the cr^^^ ^ jo lo '^bh plots 
for the CAIXAvar sample, respectively. As with the Rev sample 
plots, the excess variance is highly correlated with Mbh (with 
probability as high as 99.999 %, see Tab.|5]l. Moreover, the best 
fit lines over the entire CAIXAvar sample (solid lines in these 
panels) have slopes consistent within the errors with -1 (see 
Tab.|5]), and with the value obtained for the Rev sample. The nor- 
malisation is systematically higher in the case of the CAIXAvar 
sample. But even the largest difference of 0.4, in the case of the 
'^rms40 ^BH plots, is significant at just the ~ 2.2cr level. 



The scatter of the data around the best-fit lines in the 
CAIXAvar sample is larger than the scatter in the Rev sample. 
We found that crs(,attei-,80,40,20,io = 0.62, 0.73, 0.72 and 0.68. These 
correspond to an average scatter by a factor of ~ 5 in linear 
space. The easiest explanation is that the Mbh estimates for the 
objects in the CAIXAvar sample have an uncertainty larger than 
the uncertainty of the reverberation Mbh estimates. On the other 
hand, it is possible that the sources in the CAIXAvar and Rev 
samples do not sample exactly the same AGN population. 

To investigate this possibility, we used the cr^^^ for the 20 ks 
intervals. Fifty sources in the CAIXAvar sample are signifi- 
cantly variable on this time scale and also have Mbh estimates. 
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[OLS B] 





Fig. 4. Upper left, upper right, lower left and lower right panels show the 0"^^^,^,,, cr^^^^^, cr^^^o ^"'l lo for the 

CAIXAvar sample. The best fit relationships (see Tab.|5]l and combined 1-cr error are plotted as in Fig.|3] 



Thirty two of these objects are unique to CAIXAvar. Similarly, 
21 sources in the Rev sample show significant variations in at 
least one of their 20 ks segments. Figure |5] shows the Mbh and 
accretion rate distribution (left and right panel, respectively) of 
the 21 Rev sample sources which are variable within 20 ks (in 
blue) and of the 32 unique CAIXAvar sources, variable on the 
same time scales (light green). Application of the KS test in- 
dicates that Mbh and accretion rate distribution of the Rev and 
CAIXAvar samples are identical (a KS-test indicate that the two 
distributions are drawn from the same population at more than 
90 % probability). 

This result implies that both the CAIXAvar and Rev sam- 
ples are representative of the same AGN population. Therefore, 
the observed higher scatter of the cr^^^^ measurements around 
the best-fit models in the case of the CAIXAvar sample is pri- 
marily due to the less accurate Mbh measurements of the non- 
reverberation sources. 



6.1.3. Comparison with tine results from PSD analysis of 
the AGN light curves 

All the best fit slopes of the correlation between o-^^^^ and Mbh 
for the Rev and CAIXAvar samples are consistent (within ~ 2cr) 
with the -1 value (Tab. [3] and Tab. |5]l. This is the expected 
slope in case of an ubiquitous PSD shape with a -2 slope in 



the frequency range over which the excess variance is com- 
puted. In fact, we can predict the expected ct^^^-Mbh relation 
for the AGN in the two samples following Gonzalez-Martin et 
al. (2011). 

Panels (a) and (b) in Fig. |6] show the ct^^^^^q vs. Mbh plots 
for the Rev sample objects. Red and blue filled circles indicate 
sources with accretion rate, ih, in the range (0.0098 - 0.035) and 
(0.42 - 0.86), respectively. The average accretion rate, < ih >, of 
the objects in these two groups is 0.016 and 0.64, respectively. 
The red and blue lines in panel (a) of Fig. |6] indicate the ex- 
pected cr^^^ vs. Mbh relation according to the "case A" model 

of Gonzalez-Martin et al. (201 ifl i.e. using their equations (7) 
in the case when m = 0.016 (red line), and m - 0.64 (blue line). 
We refer to this case as "Model A" hereafter. 

Although the slope of the model lines is -1, hence they de- 
scribe well the decreasing trend of cr^^^ with increasing Mbh, 



The "case A" model of Gonzalez-Martin et al. (2011) is based on 
the results from recent detailed PSD analysis of long X-ray light curves 
of a few AGN. It assumes that the AGN PSD has a -1 slope up to a 
break frequency, Vb, , above which the slope steepens to -2. This break 
frequency decreases with increasing Mbh and increases proportionally 
with the accretion rate, as in McHardy et al. (2006), while the PSD am- 
plitude, defined as the product of the PSD value at the break frequency 
times the break frequency itself, is the same for all objects, and equal to 
0.02. 



8 



G. Pond et al.: CAIXA: a catalogue of AGN in the XMM-Newton archive III. Excess Variance Analysis 




Log(MBH) 



Fig. 5. (Left panel) The Mbh distribution of the sources which 
are variable in 20 ks intervals and belong to the Rev sample 
(blue histogram), and to the CAIXAvar sample only (light green 
histogram). (Right panel) The accretion rate distribution of the 
same sources. 




Log. of Mg^ (in Solar mass units) 



Fig. 6. (Upper left): cr^^^ 4q vs. Mbh for the Rev sample. Red and 
blue filled circles indicate sources with ih in the ranges (0.0098 

- 0.035) and (0.42 - 0.86), respectively, thus with average ac- 
cretion rates of < m >= 0.016 and 0.64, respectively. While the 
black filled circles show the cr^^^ ^ for the sources with differ- 
ent accretion rates. The red and blue lines indicate the expected 
Model A relations for the two accretion rates. (Upper right): 
CAIXAvar (20 ks) and Model A predictions (color code as be- 
fore). Lower left and right: Model B predictions. 

this model also predicts that, at a given Mbh, cTmis should in- 
crease considerably with m, contrary to what is observed. The 
blue and red points in Fig. |6] suggest that, at a given Mbh, both 
the high and low accretion rate objects have the same excess 
variance. In fact, it is because of this effect that the scatter of the 
Rev sample "cr^j^j.^,, vs. Mbh" relation is so small. 

The discrepancy between the model predictions and the 
observed cr^ .„ vs. Mbh relations becomes even more evi- 

rms,40 

dent when we consider the CAIXAvar cr^^^ 20 ^bh relation 
(shown in panels (c) and (d) of Fig. |6|. In these plots we have 
considered all the positive CAIXAvar cr^^^ 20 values (irrespec- 
tive of the magnitude of their error). Red and blue circles in these 
panels indicate the CAIXAvar objects with m in the range (0.035 

- 0.076) and (0.42 - 0.96), respectively. The red and blue lines 
in panel (c) indicate the Model A curves for m - 0.06 and 0.67 




Fig. 7. Mbh measured from the cr]^^^ vs. Mbh best fit rela- 
tion of the Rev sample, as a function of the Mbh estimated in 
other ways. Black, blue and green points show the results from 
•^rms 20 40 10' respectively. The blue line, in Fig. |2]corresponds to 
the one-to-one relation between the two Mbh estimates, and the 
red lines indicate a dispersion of a factor of 5, along the x-axis. 

(the mean accretion rate for the objects indicated with the red 
and blue circles, respectively). The model A prediction slightly 
underestimates the observed relation for the low accretion rate 
objects. A PSD amplitude larger than 0.02 (i.e. the value we have 
adopted in this work), could improve the agreement between the 
Model A predictions and the data for the low accretion rate AGN 
in CAIXAvar In this case, however, the disagreement between 
the Model A relation and the data for the high accretion rate 
AGN (blue line and blue circles, respectively) would increase 
even more. 

This discrepancy seems to suggest that one of the two as- 
sumptions in Model A is wrong. Thus, that the PSD either does 
not depend on accretion rate or its normalisation is not the same 
in all AGN. However, the PSD break frequency does depend on 
the accretion rate (McHardy et al. 2006; Koerding et al. 2007). 
This led us to investigate the possibility that the PSD ampli- 
tude is not the same in all AGN, but depends on m instead. 
A decreasing PSDamp with increasing m could counterbalance 
the expected increase of the excess variance with increasing m 
(which is mainly due to the subsequent increase of the break 
frequency). For that reason, we fitted the (c^n,s4o, Mbh) data of 
the Rev sample, using again the "case A" model of Gonzalez- 
Martin et al. (2011) but assuming that PSD^mp - AmT^. The 
best-fit results are as follows: - 46.7 for 17 degrees of free- 
dom, A = 0.003+° ™J2^ and j6 = 0.8 ± 0.15 (eiTors indicate the 
90% confidence region for a single interesting parameter). The 
red and blue lines in the bottom panels of Fig. |6] indicate the 
predicted ct^^^^q vs. Mbh relations for the Rev and CAIXAvar 
low and high accretion rate objects, using the best-fit parameter 
values we mentioned above (hereafter we refer to this case as 
"Model B"). The agreement between the Model B predictions 
and both low-m and high-m Rev sample objects has been im- 
proved, although for the high accretion rate objects a significant 
discrepancy between the model curve and the data (see the bot- 
tom left panel of Fig. |6ll is still present (this being the reason 
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Fig. 8. Upper left, upper right and lower left panels show the cr^^, g,,' ^rms,40' o"™s,20 ^rms.io 
CAIXvar sample. The best fit relationships (see Tab.|5]l and combined l-cr error are plotted as in Fig. [3] 



vs. accretion rate relation for the 



Name 


Log(MBH)opt 


Log(MBH)vai 


MRK335 


7.15 ±0.12 


6 75+" "* 



Table 4. List of AGN for which Mbh has been measured through 
variability (Log(MBH)vai) and compared to Mbh estimated in 
otherwhise. 



for the large best-fit value). However, when we consider the 
CAIXAvar cr^,-„j, -Mbh relation, the model B curves agree very 
well with both the high and low accretion rate objects (panel (d) 
in Fig.|6]). 

Furthermore, the Model B lines can also explain the fact 
that the observed cr^^^ vs. Mbh plots are well fitted by a single 
line, with the spread of the points around the best-fit line be- 
ing very small (at least in the Rev sample case). At Mbh higher 
than ~ IO^Mq, Model B predicts similar excess variances, ir- 
respective of the accretion rate of the objects. However, when 
Mbh < IO^Mq, we expect a difference up to ~ 5 in the observed 
(T^j^j. values of low and high accretion rate objects, with higher 
accretion objects showing smaller excess variance. 



6.1.4. Weighing AGN black Inoles using excess variance 
measurements 

In Section 6.1.1. we commented that the variability Mbh esti- 
mates, using the best-fit lines for the Rev sample, should be as 
accurate as the reverberation mapping ones. It seems then ap- 
propriate to use the best-fit relations for the Rev sample (with 
the parameter values as listed in Table 1) and the excess vari- 
ance measurements of the CAIXAvar sources, to estimate their 
Mbh- These estimates should have smaller uncertainty than the 
uncertainty of the present estimates which are mainly based on 
single epoch spectra. Table |4] lists the variability Mbh estimates 
for 55, out of the 161 AGN, in CAIXAvar (for 6 of these AGN 
no other Mbh estimate was found in literature). For the remain- 
ing objects, we are able to provide upper limits on their Mbh, as 
follows. 

We first measured Mbh from cr^^^^ 20 ^'^^ sources with 

at least one variable 20 ks interval (black points in Fig.|7]|. There 
are 54 such AGN in CAIXAvar (47 measurements and 7 upper 
limits). For each one of these objects, we estimate Mbh using 
their cr^_, and the best-fit cr^_, vs. Mbh relation for the Rev 
sample. The uncertainties on Mbh are estimated propagating the 
•^rms 20 uncertainty, only. Whenever the resulting uncertainty is 
smaller than 0.4 dex (estimated to be the uncertainty of the Rev 
best-fit relation itself), the latter must be used. We then followed 
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the same procedure and used the cr^ .„ and cr^ ,„, to estimate 

rms,40 rms,10' 

the Mbh for further 2 objects (one of which is an upper limit) 
and 9 more AGN (two of which are upper limits), respectively 
(blue and green points in Fig.|7]i. 

In this way, we have been able to measure the Mbh for a total 
of 55 AGN (out of 161 AGN in CAIXAvar) and provide strin- 
gent (90%) upper limits for 10 more objects. Out of these 65 
AGN, 44 have neither reverberation mapping nor stellar veloc- 
ity dispersion Mbh estimates. The remaining AGN, for which 
we detect no variability in any of the intervals analysed, are 
mainly weak sources with short exposures. However, for some 
of these sources the data are of good quality, thus indicating that 
the sources have a small variability amplitude on the time scales 
considered here. Nevertheless, even for these objects, we used 
their 90% upper limit on cr^^^^ (which is positive in at least one 
interval), and the method outlined above, to provide a lower limit 
on their Mbh- These limits are also listed in Tabled 

Figure |2] shows a plot of the variability Mbh measurements 
as a function of the Mbh estimated in other ways. Black, blue 
and green points show the results from cr^^ 40 lO' respectively. 
The blue line, in Fig. |7] corresponds to the one-to-one relation 
between the two Mbh estimates, and the red lines indicate a dis- 
persion of a factor of 5, along the x-axis. The vast majority of the 
objects lie within the red lines. This result suggests that, since 
the uncertainty associated with the current Mbh estimates for 
most of the CAIXAvar AGN are of the order of ~ 5 (see dis- 
cussion in Section 6.1.2 above), the variability estimates we list 
in this work should be considerably less uncertain, as otherwise, 
we would expect more objects to lie outside the region defined 
by the red lines in Fig. |7] 

6.2. cr^s accretion rate 

FigurelS] shows the variability (T^ gA,cr^ ^n,cr^ o„ and cr^ ,„ 

' — ' rms,80' rms,40' rms,2() rms,10 

VS. accretion rate relation for the CAIXAvar sample. The best-fit 
parameters for each plot in Fig. [8] are listed in Table |5j together 
with the probability of significant correlations. 

The plots in Fig. [8] suggest the presence of a trend where 
higher accretion rate sources appear to be more variable. But 
there is also a significant scatter around the best-fit relation 
(shown with a solid line in the same plots). Not surprisingly, 
given the large scatter of the points around the best-fit lines, the 
probability of a significant correlation between excess variance 
and accretion rate is small (even smaller than 90% in almost all 
cases). Therefore, strictly speaking, our results suggest that there 
is no significant correlation between the variability amplitude 
and the accretion rate in CAIXAvar 

The scatter of the points around the best-fit lines in Figure[8] 
may be caused by the strong variability vs. Mbh correlation in 
the CAIXAvar sample. If less massive BH are intrinsically more 
variable than higher mass BH, they should populate the higher 
part of the variability vs. accretion rate plot at any given accre- 
tion rate. Thus, even if there exists a variability vs. accretion rate 
correlation, we should expect a significant scatter in the plots 
shown in Fig. [8] simply because of the large Mbh range of the 
objects in our sample. 

In order to eliminate the Mbh dependence of 0-^^^,^, we con- 
sidered the product cTr^s 20 ^^bh- Given the fact that the best- 
fit slope of the cr^^^o vs. Mbh relation is consistent with -1 
(see Table |5]i, this product should result in a quantity inde- 
pendent of Mbh- The bottom panel in Fig. |9] shows a plot 
of the cr^s 20^^BH vs. accretion rate relation for all AGN in 
CAIXAvar (irrespective of the amplitude of their excess variance 




Bol' '-Edd 



Fig. 9. The ct^^sIO^^^bh 

vs. ih relation for the CAIXAvar (20 
ks) sample. Red and blue circles indicate sources with Mbh in 
the range 1.7-4.6x1 Q^Mq and 3-8 . 9 x 1 O'^Mq . Red and blue lines 
indicate the expected relations for mean Mbh of 2.5 x 10^ and 
5 X IO^Mq, respectively. Dashed and solid lines show the Model 
A and Model B predictions. 



uncertainty). Clearly, a large scatter in this plot is introduced by 
the uncertainties on the Mbh, and of the bolometric luminosity 
estimates for each object (as before, LboI is estimated using the 
prescription of Marconi et al. 2004), which should increase the 
uncertainty of the LBoi/LEdd values. In any case though, and in 
agreement with previous studies (O'Neill et al. 2005; Gierlinski 
et al. 2008; Nikolajuk et al. 2009; Zhou et al. 2010), no correla- 
tion between ct^^^^xMsh and accretion rate appears. We therefore 
conclude that cr^^^^ does not appear to correlate with accretion 
rate. 

To investigate this issue further, the red and blue circles 
in both panels of Fig. |9] indicate the cr^^^ vs. m (top panel) 
and cr^j^j, xMbh vs. m data (bottom panel) for AGN in the 
CAIXAvar sample with Mbh in the range 1.7-4.6x10*' Mq and 
3-8.9x10^ Mq, respectively. The solid (dashed) red and blue 
lines indicate the Model B (Model A) curves for an AGN with a 
Mbh of 2.5 x 10* and 5 x IO^Mq (these values are equal to the 
mean Mbh of the points indicated with the red and blue circles 
in Fig.|6]l. 

The top panel in Fig. |9] shows that, because of the cr^^^ vs. 
Mbh dependence, neither Model A nor Model B predict a univer- 
sal cr^i^j. - m relation for AGN. We do expect a significant scatter 
in the cr^^^^ vs. ih plot, because, at any given accretion rate, the 
excess variance should be different for AGN with different Mbh- 
In fact, it should increase strongly with decreasing Mbh, due to 
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Fig. 10. The cr^^_, vs. bolometric luminosity, and the <t^^^ 20 ^ '^bh vs. bolometric luminosity plots for the CAIXAvar sample (left 
and right panels, respectively). Lboi,44 indicates the bolometric luminosity in units of 10'^'* erg s 



the "PSD break time scale - Mbh" relation of McHardy et al 
(2006). On the other hand. Model A predicts a cr^^^ oc m rela- 
tion, while Model B predicts a much shallower relation of the 
form cr^^j, oc m*^^, due to the decrease of the PSD amplitude with 
increasing accretion rate. The McHardy et al (2006) relations 
predict that the PSD break-time scale becomes smaller than 20 
ks (thus entering in the frequency window in which the cr^^^ is 
computed) for an AGN with Mbh = 2.5 x 10'' Mq (red curves in 
the top panel of Fig. |9|l when rii ~ 0.05. As a result, the Model 

A, with constant PSD^mp, predicts that the excess variance of 
the small Mbh should increase linearly with m until m = 0.05, 
then bend until the PSD break time scale becomes even smaller 
than 250 s. At that point the excess variance should flatten, be- 
coming constant. While the more massive AGN should always 
have the PSD break time scale longer than 20 ks, thus no flat- 
tening is expected. When we consider the predictions of Model 

B, on top of this trend with accretion rate, we also have to con- 
sider the anti-correlation between PSD^mp and ih. In Model B, 
the PSD normalisation vs. accretion rate relation will reduce 
the expected trends predicted by the Model A relation by 0.8 
(PSDatnp oc m""^), producing, i.e. for the small Mbh, a flatter 
relation at low m and a decreasing variability at high m. 

The main effect of multiplying the excess variance by the 
Mbh is to eliminate of the Mbh dependence in the cr^|^j.vs. ih 
relation. In fact, as expected, the predicted values of ct^^^xMbh, 
for different Mbh do overlap both in Model A and B scenarios, 
at least when the PSD break time scale is not entering in the 
frequency window in which cr^^^^ is computed (i.e. lower than 20 
ks). 

The Model A (dashed curves) predictions in the bottom 
panel of Fig. |9] clearly fail to reproduce the behaviour of both 
the low and high ih AGN. The low ih sources show cTj^j^j, xMbh 
values which are significantly larger than the Model A curve. At 
high ih, although the model correctly predicts a flattening of the 
relation for the small Mbh, it systematically overestimates the 
cTj^sXMbh values. As before, a PSD^mp lower than 0.02 could re- 
sult in model curves which are consistent with high m small Mbh 
AGN, however, in this case, the disagreement between Model A 
curve and the low ih AGN (red dashed line) will get even worse. 
On the other hand, the data for all the low ih objects are now 
consistent (within 2cr) with the Model B curve (solid line), i.e. 
the solid blue line shown in the bottom panel of Fig. |9l The solid 



red line (in the same panel) is also in agreement with the data for 
the low Mbh objects. 

Therefore, the apparent lack of correlation between cr^^^ 
(when "normalized" to Mbh) and ih strongly supports the hy- 
pothesis that the PSD amplitude in AGN decreases with increas- 
ing accretion rate. We conclude that the cr^^^vs. m relation for 
CAIXAvar is consistent with the McHardy et al. (2006) Vbr m 
correlation, but only if the PSD normalisation decreases with in- 
creasing accretion rate. 



6.3. cr^njs vs. luminosity 

The left panel of Fig.[TO]shows the plot of cr^^^ 20 '^boi for the 
CAIXAvar sample. Table|5]lists the best linear fit (in the log-log 
space) in the case of the cr^^^ 20 10 '^boi relations, as well 
as the probability of a significant correlation between these two 
quantities. Variability and luminosity are well correlated with 
correlation probabilities higher than 99.9 %.The best-fit slopes 
are also consistent with -1, in all cases. 

The right panel in Fig. [TOl shows the o-^i^s 20^^^^ ^boI 
relation. Contrary to the data plotted in the left panel of the 
same figure, no significant correlation is observed any more be- 
tween cr^^j, xMbh and luminosity (see also Tab. |5]l. This result 
suggests that the apparent correlation between cr^^^ and LboI is 
driven mainly by the cr^^^ vs. Mbh correlation. This result is also 
consistent with the suggestion that the variability vs. luminosity 
relation in AGN is a byproduct of the variability vs. Mbh relation 
in the same objects (see e.g. Papadakis 2004). 

To investigate this issue further, we computed the expected 
Model A and Model B cr^,^^ vs. LboI relation for various accretion 
rates. The red and blue lines in the top and bottom panels of Fig. 
[TT]indicate the Model A (top panel) and Model B (bottom panel) 
relation for AGN with m = 0.06 and m - 0.67, respectively. 
Black circles in both panels indicate all the cr^^^ vs. LboI data 
for the CAIXAvar sample, and red/blue circles indicate the data 
for sources with an accretion rate between (0.0098 - 0.035) and 
(0.42 - 0.86), respectively (as we did in Section 6.3.1). In effect, 
the model lines in this figure are identical to the model lines 
plotted in panels (c) and (d) of Fig.|6] except that, for each Mbh, 
we have estimated Lem, and hence LboI, through the relation: 
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LboI = tiihEdd (note that, due to the multiplication of LEdd with 
m, even the Model B curves are well separated in this case). 

Although the agreement between Model A curve prediction 
with the low ih objects is rather good, this is not the case with the 
model prediction for the high accretion rate objects (all points 
lie below the blue line in the top panel of Fig. fTTT i. On the other 
hand, the Model B curves agree very well with both the low and 
the high m data points. We therefore conclude that the apparent 
variability vs. luminosity relation in the CAlXAvar sample is 
a byproduct of the variability vs. Mbh relation, and it can be 
explained well if the PSDamp decreases with accretion rate as we 
discussed in Section 6.3.1. 

6.4. cr2„, VS. Hp 

The left panel of Fig.[T2]shows the cr^^^ 20 ^^^^ Width 

at Half Maximum (FWHM) relation. As observed by McHardy 
et al. (2006), a strong anti-correlation is present in CAlXAvar 
(the best-fit results, as well as the probability of a correlation 
between variability amplitude and FWHM are listed in Table 
|5]|. However, it is well know that Hp FWHM strongly depends on 
Mbh- In fact, Bianchi et al. 2009 show a very strong correlation 
between Mbh and H^ for the objects in the original CAIXA sam- 
ple. This is not surprising, in fact the width of the Hp line is often 
used to estimate Mbh in AGN. Therefore, the anti-correlation be- 
tween variability and Hp FWHM that we observe, is almost cer- 
tainly a byproduct of the variability vs. Mbh relation. The right 
panel of Fig. [12] shows that the variability vs. width of the H^ 
line correlation disappears when the dependance on the Mbh on 
excess variance is taken off by multiplying crl^^ with the Mbh- 
This is true also for the excess variance measurements from the 
80, and 40 ks segments (see Table|5]|. 

6.5. crl^^ VS. spectral index T 

The left panel of Fig.[T3]shows the cr^^^ 20 2-10 keV spec- 

tral index relation for the CAlXAvar sample (see Bianchi et al. 
2009 for a description of how the spectral index is calculated). 
The variability is strongly correlated with F. This is also the case 
for cP- „„, cr^ .„, cP- -„ and the cr- , „ vs. F relations (see Tab. 

rms,ftO' rms,40' rms.zu rms,10 ^ 

|5]). In fact, the probability that these correlations are significant 
are all above 98.5 % and reaches values larger than 99.99 %. 

This correlation (see also Turner et al. 1999) is rather surpris- 
ing and not easy to understand. As before, we also investigated 
the cP^^ 20 io^^Mbh vs. F relations. When the variability am- 
plitude is "normalized" to Mbh, the correlation disappears (see 
right panel of Fig. [13] and Tab. O. However, even this result is 
not easy to understand, as there are no "a priory" direct physical 
connection between the Mbh and the 2-10 keV spectral index. 
In fact, although some authors claim the presence of a correla- 
tion between spectral index and Mbh, only a marginally signifi- 
cant correlation between these quantities is observed in CAIXA 
(Bianchi et al. 2009). To be scrupulous, we repeated the study of 
this correlation in the case of the CAlXAvar sample (due to the 
slightly different definition of the CAlXAvar sample compared 
to CAIXA). Again, only a marginal correlation is present. 

Several authors have argued in the past that the 2-10 keV 
spectral index is not correlated with the Mbh, but with the ac- 
cretion rate (see e.g. Porquet et al. 2004; Piconcelli et al. 2005; 
Shemmer et al. 2006; Saez et al. 2008; Papadakis et al. 2009; 
Sobolewska & Papadakis 2009; Risaliti et al. 2009; Scott et 
al. 2011). This argument is strengthened by the fact that the 
same spectral index vs. accretion rate correlation is observed in 
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Fig. 11. Top panel: The cPj^^20 ^BoI relation for the 
CAlXAvar (20 ks) sample. Red and blue circles indicate sources 
with mean m = 0.06 and 0.67, respectively. The solid lines show 
the Model A predictions for different accretion rates. Bottom 
panel: Model B predictions. 



BHB as well (see e.g. Wu & Gu 2008). If this is ti'ue, then we 
can explain the crj^^^ vs. F relation as follows: Sobolewska & 
Papadakis (2009) have argued that F oc m" ', while we observe 
that cr^ oc m (see Table |5]l. Therefore, we would expect to ob- 
serve cr^j^j, oc F'", perfectly consistent with our observations (see 
Table|5]l. 

In Fig. [14] we plot again the cr^^^ vs. F data for all the 
CAlXAvar sources (irrespective of the magnitude of their er- 
ror), in the log-log space. Red and blue circles indicate the 
data for AGN with 1.7 x 10^ Mq < Mbh < 4.6 x 10^ Mq 
and 3 X 10** Mq < Mbh < 8.9 x 10** Mq, respectively, ex- 
actly like in Fig. [6] The red and blue lines indicate the Model 
B cr^j^, - F predictions for an AGN with Mbh = 2.5 x lO*" and 
5 X 10** M0 (i.e. the mean Mbh of the points indicated with the 
red and blue circles). For each of the two Mbh values we have es- 
timated the excess variance for various accretion rates assuming 
that PSDamp oc m as explained in Section 6.1.3. Then, instead 
of plotting cr^ms as a function of m (as we did in Fig. |9]l, we plot 
it as function of F, using the Sobolewska & Papadakis (2009) 
relation mentioned above. The agreement between the data and 
the model predictions is rather good. This result suggests that 
the cr^ms ~ r relation we observe (Fig.[T3Ti is due to the o-^^^ vs. 
ih and F vs. ih relations. When "normalized" to Mbh, the excess 
variance is no more clearly correlated with F, because, according 
the Model B predictions, it is expected to be weak (with a slope 
of just ~ 0.2), and with a considerable scatter at high accretion 
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Fig. 12. (Left panel) cr^^^^o vs. FWHM H^. A strong correlation is present. (Right panel) cr^j^^2o ^ Mbh vs. FWHM H^. The 
correlation disappears. Hy62ooo indicates Hy6 in units of 2000 km s ^ 
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Fig. 13. The cr^^jn ^nd cr^^jo ^ ^bh vs. 2-10 keV spectral index correlations (left and right panel). Fu indicates F in units of 1.7. 



rates, where the excess variance may increase with increasing m 
for AGN with large Mbh, but it may also decrease with increas- 
ing ih for lower Mbh AGN. 

7. Are NLS1 intrinsically more variable than other 
AGN? 

Several works in the past have shown that NLS 1 appear to have 
higher variability than normally expected for their Mbh (see e.g. 
Leighly 1999; Turner et al. 1999; Nikolajuk et al. 2004; 2009). 
To investigate this issue, we compared the excess variance of the 
NLSl and BLSl in CAlXAvar. The upper left panel of Fig.fTSl 
shows the distribution of Log(cr^j^^ ^q) foi" the NLS 1 (light green) 
and the other AGN (BLSl, in blue) in the CAlXAvar sample. 
Clearly, the NLS 1 appear to be more variable than BLS 1 . The 
mean value of the distributions are < Log(cr^^^ 2())> = -1-69 and 
< Log(cr2^j 2o)> = -2.70 for the NLSl and BLSl, respectively. 
Application of the Student's t-test indicates that they are signifi- 
cantly different (null hypothesis probability, NHP= 3 x 10"^'). 
At the same time, the upper middle and upper right panels 
of Fig. [13] shows that NLSl have significantly smaller Mbh 



( <Log(MBH,NLSi) >= 6.86 and <Log(MBH,BLSi) >= 7.85; 
NHP= 2 X 10"^). They also appear to have higher accretion rates 
(upper right panel in Fig. [TSj, although the difference between 
NLSl and BLSl ih is less significant ( <Log(mNLSi) -0.47 
and <Log(mBLSi) >= -0.89; NHP= 2.2x10"^), just a ~ 1-sigma 
effect. 

The lower left panel of Fig. [15] shows the cr^^^^ 20^'^BH dis- 
tribution for the NLS 1 and BLS 1 as well. They look very sim- 
ilar and in fact the application of the t-test indicates no signifi- 
cant difference between the mean of the two distributions (NHP 
> 0.05). We also test that, given the number of points in each 
sample, we would be able to detect with high significance (more 
than 3cr) a difference in the mean of the two distributions of 
the order of 4-5. This indicates that, once the Mbh dependence 
is taken into account, NLS 1 have variability amplitudes, to first 
approximation indistinguishable from the ones of BLS 1 . Thus 
the difference in the variability properties of NLS and BLS 1 is 
mainly due to the fact that NLS 1 host a smaller Mbh, on average, 
when compared to BLS 1 . 

In order to investigate whether the accretion rate also con- 
tributes to the difference in the variability amplitudes of NLS 1 
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Fig. 14. The cr^^^^o vs. 2-10 keV spectral index correlation. 
Solid line show Model B predictions for different Mbh- 

and BLS 1 we selected a subsample of NLS 1 and BLS 1 which 
have the same Mbh, on average. The upper middle panel of Fig. 
[TSl indicates that several NLS 1 and BLS 1 have Log(MBH) be- 
tween 7 and 7.6. We note that MRKllO belongs to this Mbh 
range. The classification of this AGN is quite uncertain. In fact, 
according to the FWHM of Hp = 1760 km s it would be clas- 
sified as a NLSl, however it lacks all the other features typical 
of NLSl. For this reason, detailed studies of its optical spec- 
trum suggest that MRKllO is a BLSl with relatively "narrow" 
broad lines (Veron-Cetty et al. 2007). Thus, we decide to exclude 
MRKl 10 from the following analysis. Both the NLS 1 and BLS 1 
distributions in this narrow Mbh range have mean Log(MBH) of 
7.37. 

The lower middle panel of Fig. [15] shows the excess vari- 
ance distribution of the NLS 1 and BLS 1 within the selected Mbh 
bin mentioned above. Despite the small number of objects con- 
sidered here, the distributions suggest that NLSl show higher 
variability (KS and t-test more than 98.7 % probability; NHP 
= 1.3 X 10"^; the mean excess variance being < Log((r^^^ jo)'* - 
-2.04 and -2.70 for NLSl and BLSl, respectively). This differ- 
ence cannot be ascribed to the Mbh, thus either the NLS 1 are dif- 
ferent from the other AGN or the difference is due to their higher 
accretion rate (< m >nlsi= -0.65 and < m >blsi= -1.03, al- 
though not significantly different, NHP> 0.05). Thus, to test if 
the enhanced variability is due to the accretion rate dependence, 
we multiply the excess variance for Mbh and divide for the ac- 
cretion rate to the expected power (m' as in Model A and m*'^ as 
in Model B). The right lower panel of Fig. [IS] shows that, once 
that the accretion rate is taken into account (showing the case B 
of m"^), the distribution and the mean of the variability of NLS 1 
and BLSl are perfectly consistent (NHP > 0.05 in both cases). 

8. Discussion 

We have studied the X-ray variability properties of a large (161 
AGN) sample of radio quiet. X-ray unobscured, bright AGN. 
We estimated their excess variance on various time scales (80, 
40, 20 and 10 ks) and various energy bands (0.3-0.7, 0.7-2 and 
2-10 keV), using XMM-Newton data and taking full advantage of 
the high sensitivity, long orbit and wide energy range of XMM- 
Newton. We have systematically investigated the relation be- 
tween the variability amplitude of these objects (as parametrized 
by cr^^J and other source parameters like, Mbh, accretion rate. 



optical emission line width, and X-ray spectral slope. The AGN 
of the CAlXAvar sample, with well measured cr^,^^, mainly be- 
long to the local universe (z< 0.2, with just 3 objects at interme- 
diate redshift, see Fig. 1) and have relatively modest bolometric 
luminosity (LboI ~ 10'""'*^ erg s"'). Thus, the observed relations 
and variability properties are, strictly speaking, characteristic of 
the AGN of the local Universe. A larger sample of higher redshift 
and luminosity AGN is needed in order to confidently extrapo- 
late the AGN variability properties observed in the local universe 
to higher redshift. We discuss below some of the implications of 
our results. 



8. 1. cr^^s in various energy bands 

One of the main results of our work is that the variability ampli- 
tude is comparable in the hard, medium and soft energy bands, 
at least on time scales < 40 ks. Typically, different spectral com- 
ponents contribute to the AGN X-ray spectra at different energy 
bands. At energies above 2 keV, the spectrum is generally dom- 
inated by the continuum power-law component. At lower ener- 
gies a "soft excess" component appears in many objects, on top 
of the extrapolation of the power law to lower energies. If this 
component is constant (or less variable than the continuum), then 
we should expect cr^ r, < cr^ t mi, v, but this is not the case. 

rms.soft rms,2-10keV' 

Furthermore an absorbing material (either neutral or ionised) is 
also present in many radio quiet AGN. This material acts as a fil- 
ter, absorbing mainly the soft and medium energy band photons, 
and may vary in ionisation, responding to the continuum varia- 
tions, and/or in its column and covering fraction of the source. 
If these variations were to happen on time scales less than ~ 
half a day then we would expect to measure a larger variability 
amplitude in the medium and soft energy bands, since they are 
more affected by absorption, compared to the hard band. This 
is contrary to what we observe. Finally, if there were signifi- 
cant, large amplitude intrinsic spectral slope variations with a 
pivot point at energies lower than 0.3 keV, then we would ex- 
pect cr^ T ,ra, > cr^ ft J , and the opposite trend if 

^ rms.z-lukeV rms.solt, or medium' 

the pivot point were at energies higher than 10 keV. Again, this 
is contrary to our results. Consequently, our results strongly in- 
dicate that, on time-scales less than 40 ks, intrinsic spectral slope 
variations or absorption variations, must be of small amplitude 
so that they do not contribute significantly to the observed ^m.x: 
variability of the sources in the different energy bands. Of course 
we know that the rms-spectra of AGN are not flat (e.g. Ponti et 
al. 2004; 2006; 2007; 2010; Gallo et al. 2004; 2007; Gierliriski & 
Done 2006; Goosman et al. 2006; Petrucci et al. 2007; Larsson 
et al. 2008; Zoghbi et al. 2010). In other words, the components 
mentioned above, do affect the observed variability at different 
energy bands to some extent. However, our results show that, 
to first order, the main driver for the observed^MJc variations in 
AGN at all bands, between 0.3-10 keV, is the continuum nor- 
malisation variability. 

On 80 ks time-scales a deviation, although not more signifi- 
cant than ~ 2 cr, from the 1-to-l correlation appears. In particu- 
lar, the more variable (i.e. the smaller Mbh) AGN appear to have 
a higher soft and medium variability, when compared to the hard 
band. This trend had already been noticed in the past as well (e.g. 
Nandra et al. 1997; Leighly 1999; Markowitz & Edelson 2001). 
If real, this result is most probably due to a complex variation 
of the PSD normalisation and slope with energy (McHardy et al. 
2004; Utdey & McHardy 2005). Alternatively, this trend might 
indicate some sort of triggering of the soft/medium band vari- 
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ability occurring only on time-scales longer than a given time- 
scale and-or characteristic length. 

If this deviation from the 1-to-l correlation in the cr^ „,i 

rms,oU 

plots is real, then we would expect it to appear even stronger at 
longer time-scales. Indeed, Markowitz & Edelson (2001) mea- 
sured the cr^^j on time-scales of 300 days (~ 26 Ms) with time 
bins of roughly 5-13 days (~ 800 ks) for 9 Seyfert 1 galaxies. 
On such long timescales, the authors find that all the sources ex- 
hibit stronger variability in the softer X-ray (the 2-4 keV band, 
in their case) than in the hard one (7-10 keV band). When the 
o"Ls,7-io keV is plotted vs. the o-^^, 2-4 keV a highly significant de- 
viation from the 1-to-l correlation is observed (best fit slope 
being 0.724 + 0.041; Markowitz & Edelson 2001). This result 
strongly supports the idea of a triggering of the spectral variabil- 
ity at a timescale of about 80 ks (at least for the low Mbh AGN). 
For this reason we discuss below some ideas which could ex- 
plain such a behaviour. 

Absorbers change their ionisation balance (i.e. change their 
opacity at a given energy), responding to the illuminating 
sources, only on time-scales longer than their recombination 
time (on shorter time-scales the absorber has not enough time 
to react to the variation of the illuminating source). The recom- 
bination time depends critically on the density of the absorber, 
being shorter for denser clouds. Interestingly accretion disc the- 
ory suggests that smaller mass BH should have denser discs and 
environment. If this is the case it would indicate that the ab- 
sorber's recombination time scale is of the order of 40-80 ks for 
small Mbh AGN (and longer for larger Mbh AGN). A time scale 
of 80 ks corresponds to about 1 light day (roughly the distance of 
the BLR in a small mass BH) or a light crossing time of 2 x lO'* 
rg, for a BH of 10'' M© (r^ = ^). Such radius is compara- 
ble to the Compton radius where a disc wind might be formed 
(Begelman et al. 1983a,b; Woods et al. 1996; Ponti et al. 2012). 
Assuming that the excess of soft and medium variability is due 
to absorption located either in the BLR (see also Risaliti et al. 
2011) or at the Compton radius, we would naively then expect 
the absorber to be variable on time-scales larger than about sev- 
eral IO^xMbhj s (see Ki-ongold et al. 2005; 2007; Kaastra et al. 
201 1; for more accurate computations). 

Alternatively the spectral variability might be intrinsic to the 
source, with the spectral slope being correlated with flux, having 
a pivot point above 10 keV, thus producing higher soft variabil- 
ity. Such a phenomenon has been envisaged by Comptonisation 
models (Haardt et al. 1991; 1993; 1994; 1997) in which varia- 
tions of the coronal optical depth, size and-or temperature drive 
spectral changes, with the spectral index steepening with flux. 
If the time scale for significant spectral slope variations corre- 
sponds to the thermal and-or sound crossing time of a standard 
accretion disc at a few r^, this could be of the right order of 
~ 40 - 100 ks for small Mbh AGN (where the thermal and sound 
crossing time should be ~ 30 - 150 ks for an alpha disc with 
a = 0.01, H/R=0.01 at 7 rg for Mbh = 10'' Mq; Treves et al. 
1988; Edelson & Nandra 1999). 

8.2. cr^^, vs. /Wbh, bolometric luminosity and Hp 

The CAIXAvar sample shows statistically significant correla- 
tions between cr^jj^j, and Mbh, bolometric luminosity and H/S 
(McHardy et al. 2006). We believe that the latter correlations 
are a by-products of the cr\-^-^^ vs Mbh relation, which is probably 
the fundamental correlation in AGN (Papadakis 2004). There are 
two reasons why this should be the case: a) the scatter of the cr^^^ 
vs. Mbh is much smaller than the correlation of the cr^^^ with 



LboI or with the FWHM of Hj6; b) standard accretion disc theory 
(Shakura & Sunyaev 1973) predicts that all the time-scales in an 
ff-disc (like the light crossing, orbiting, thermal, sound crossing, 
radial drift, and viscous time scale), which should affect the vari- 
ability amplitude of an object, scale linearly with Mbh (Molendi 
et al. 1992; Treves et al. 1988). 

Our best-fit slope and normalisation of the ct^^^q vs. Mbh 
relation are in agreement with the results of Zhou et al. (2010), 
who considered a sample of AGN with Mbh from reverberation 
mapping only. Their best-fit slope of a^Q - -1.05 + 0.08 for 
the cr^^j, -Mbh relation is in agreement with our best-fit slopes 
of both the CAIXAvar and Rev samples. O'Neifl et al. (2005) 
computed the cr^^^^ for a large ASCA AGN sample. They also 
found a significant correlation between normalized excess vari- 
ance and Mbh, and estimated a slope of ci'40 = -0.57, steepening 
to Qr4o = -0.84 (roughly in agreement with the slope obtained 
with CAIXAvar) once a broken power law was considered. 

One of our major results is that the scatter in the cr^^^ vs. Mbh 
relation of the Rev sample is very small. It corresponds to a scat- 
ter of 2.7-3 in linear space, which is comparable to the scatter 
of 2.6-2.9 that Peterson et al. (2004) measured around the AGN 
Mbh vs. ct, relation. This result implies that there may be no 
intrinsic scatter in the cr^^^ vs. Mbh relation, and the scatter we 
observe may be entirely due to the uncertainty on the Mbh values 
from reverberation. Zhou et al. (2010) used a Nuker's linear re- 
gression method to estimate the intrinsic scatter of this relation. 
They fitted the relation taking into account the uncertainties as- 
sociated with both quantities {cr^^^ and Mbh), then they added a 
component (that would measure the intrinsic scatter) to the cr]^^ 
uncertainties until reaching ilx^ - 1- 1" this way they estimated 
an intrinsic scatter of 0.2 dex (in agreement with the Rev sample 
results that show that intrinsic scatter on cr^^^ < 0.4 dex). 

A scatter around the best fit relation of less than a factor 
of 3 is challenging for the current variability models. We know 
that the reverberation sample spans a wide range of accretion 
rates from below 1 % up to the Eddington limit. If the break fre- 
quency scales linearly with accretion rate (McHardy et al. 2006; 
Koerding et al. 2007), and the PSD amplitude is constant in all 
AGN, we would expect a scatter of more than ~ 10-20 for the 
20 and 40 ks relations at Mbh > IO^Mq (see Fig. |6]l. We found 
that the tight correlation between variability and Mbh can be ex- 
plained if the PSD normalisation decreases with increasing ac- 
cretion rate as PSDamp °^ m""-^. In this case, the combined varia- 
tions of the Vbr and PSD normalisation with accretion rate results 
in the PSD integral above the longest frequency sampled in the 
light curves (i.e. cr^„ J being only weakly dependent on accretion 
rate. Thus, the expected cr^^^ at a given Mbh is almost the same 
even for objects with a large difference in accretion rates. We 
therefore conclude that, the presence of the cr^^^^ vs. Mbh corre- 
lation, its best-fit slope and normalisation, and its small scatter, 
are in agreement with the idea that every object in CAIXAvar 
has the same PSD with low frequency slope of - 1 breaking to 
a slope of -2 at a break frequency that decreases with Mbh and 
increases with m (McHardy et al. 2006; Koerding et al. 2007), 
but only if the PSD normalization also decreases with increasing 
m. This is a new result, which has never been observed before. 
Investigation of the published PSD results indicates that this is 
indeed the case (Papadakis et al., in preparation). 

8.3. Tlie X-ray variability Mbh estimation in AGN 

The use of the cr^^^ as a tool to measure the Mbh has already 
been proposed in the past. Nikolajuk et al. (2004) provided the 
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excess variance ' MBH excess variance 7 0<Log(MBH)<7 S (excess variance * IVIBH) / IVldct^^ 



Fig. 15. (Upper left panel) The light green histogram shows the excess variance distribution of the NLSl which are variable in 20 
ks segments and present in CAIXAvar The blue histogram shows the excess variance distribution for the other (broad line Sayfert 
1) sources in CAlXAvar. The NLSl are significantly (KS test) more variable than BLSl. (Upper middle panel) Mbh distribution 
(colour code as before). The NLSl have smaller BH masses than BLSl. (Upper right panel) Accretion rate distribution. The NLSl 
are accreting at significantly higher accretion rates than BLSl. (Lower left panel) Crn,s2o^^BH distribution. Once that the Mbh 
dependence is taken into account the NLSl are not anymore distinguishable from BLSl. (Lower middle panel) Excess variance 
distribution of the AGN with 7.0<Log(MBH)<7.6. At the same Mbh, NLSl are more variable than BLSl objects. (Lower right 
panel) cr^^^^ jqXMbh/w"'^ distribution for the AGN with 7.0<Log(MBH)<7.6. The higher vai-iability of the NLS 1 objects can be due 
to their higher accretion rate. 



first recipe to measure AGN Mbh from crl^^ measurements. 
Their recipe could work for both BLSl and NLSl, but only if 
a shift of a factor of 20 in Mbh was used in the case of NLS 1 . 
With the introduction of the correction factor, the authors were 
in effect assuming a dependence on accretion rate (NLS 1 have 
in general higher accretion rates than BL AGN, see i.e. ii6.5l 
and !j7|, but rather than a continuous scaling with m (such as 
i.e. McHardy et al. 2006; Koerding et al. 2007) they were as- 
suming a bimodal dependence in the form of a correction fac- 
tor for NLSl. A few years later, Gierlinski et al. (2008) showed 
that the PSD integral above Vbi remains approximately constant 
for a given AGN and BHB (in the hard state), regardless of the 
source luminosity (thus accretion rate), unlike Vbr which corre- 
lates strongly with the latter. They argued that it is for this reason 
that the Mbh estimates from crl^^ measurements of light curve 
segments which are shorter than 1 /vbr are accurate. To further in- 
vestigate the difference between NLS 1 and BL AGN, Nikolajuk 
et al. (2009) computed the cr\-^^ for a large sample (21) of NLSl 
and confirmed that, if the recipe to estimate Mbh from cr^^, was 
deduced from the PSD scaling with Mbh only (without contain- 
ing a term related to m), a shift of a factor of 20 had to be used 
for NLS 1 . However, the authors also showed that this factor is 
consistent with being just the product of a weak dependence of 
cr^s with m. Interestingly, the measured dependence of the vari- 
ability with m was weaker (with a slope flatter a - 0.79) than the 
one expected assuming the McHardy et al. (2006) relation (with 
a - 1.4). This finding further supports the idea that the PSDamp 
lowers with m, as proposed in this work. 



Given the small scatter in the cr\-^^ vs. Mbh relation for the 
Rev sample objects, we argue that the best-fit results that we 
have listed in Tables [3] can be used to measure the Mbh for an 
AGN if any of the o"^,^^ 80 40 20 10 ^"ilu^s is available. If the cr^^^ is 
measured with an uncertainty comparable to the uncertainty of 
cr^jns for the AGN in the Rev sample, this method can potentially 
yield Mbh estimates as accurate as the ones from the reverber- 
ation mapping. In fact, we followed this approach in this work, 
and we provide new Mbh estimates for 55 AGN in CAlXAvar, 
of which 6 had no previous Mbh estimates. 47 (1 and 7) of our 
estimates were based on cr^^^ measurements from 20 ks (40 and 
10 ks) long light curves. These estimates should be on average 
as accurate as the Mbh estimates from reverberation mapping, 
and most probably more accurate than the Mbh estimated from 
single epoch spectra. Moreover, we find that this recipe can be 
applied to both BLS 1 and NLS 1 . 

Since it is easier to estimate the cr^^^ than to perform a se- 
ries of reverberation mapping measurements for an AGN, X-ray 
variability based methods should be useful in determining Mbh 
of many objects. However, we stress that one has to be careful 
when using the crlj^^ vs. Mbh best fit results listed in Table |3] 
(and in fact, any X-ray variability based method to estimate the 
Mbh in an AGN). Figure [16] shows the expected cr^^^ vs. Mbh 
relations according to our Model B prescription, for AGN with 
an accretion rate of 0.01 (red lines) and 1 (blue lines) and for 
light curves which are 80, 40, 20 and 10 ks long (from the left 
to the right panel), binned at 250 s. A shorter bin size does not 
affect the expected relation significantly, except for the smallest 
Mbh and highest accretion rate sources. To illustrate this point. 
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the dashed lines in the 3rd panel of Fig.[T6]indicate the predicted 
cr^i^j vs. Mbh relation for a light curve with a bin size of 50 s (as 
opposed to the solid lines which corresponds to a 20 ks long, 250 
s binned light curve). The relations are identical in the case of the 
low accretion rate objects (hence it is not possible to separate the 
dashed and solid red lines in this panel). On the other hand, a 
shorter bin size results in a different variability vs. Mbh relation 
for the high accretion rate objects, but only when Mbh ^ 5 x 10^ 
(see blue solid and dashed lines in the same panel). Even in this 
case though, the difference between the two relations is really 
small. 

The important thing to notice in all panels in Fig. [16] is 
that Model B predicts a small scatter in the variability vs. 
Mbh relation but only above a certain Mbh which is equal to 
3 X 1Q\ 1.8 X 10^ 10^ and 5 x lO'' Mq, in the case of the 80, 40, 
20 and 10 ks long light curves. This effect is equivalent to what 
the previous studies have argued: the Mbh estimates which are 
based on the use of cr^^^ measurements are "valid" as long as the 
length light curve used to measure the cr^^^ is shorter than 1/vbr. 
For example, if we use an <t^^^ estimate from a 80 ks long light 
curve for an object which has a Mbh smaller than ~ 3 x 10^ Mq, 
and an accretion rate larger than ih ~ 0.6 - 0.7, the break time 
scale will be shorter than 80 ks, hence, due to the decrease of 
PSDamp with increasing accretion rate, the expected excess vari- 
ance will be smaller than the excess variance of objects which 
accrete at a smaller rate. In this case, the use of the best fit re- 
sults we list in Table[3](which indicate no break of the variability 
vs. Mbh relation at low Mbh objects) may result in an inaccu- 
rate Mbh estimate. One should be careful if the use of the Rev 
sample best fit results yields Mbh estimates significantly smaller 
than the masses mentioned above. In this case, an accurate esti- 
mate of Mbh requires knowledge of the LBoi/LEdd ratio for the 
object, which may not be possible (even if LboI is known, Ledd 
requires an apriori knowledge of Mbh)- Obviously, it is prefer- 
able to use 10 ks light curves to estimate the Mbh, however, such 
a short light curve may not result in a positive excess variance 
measurement, at least for large Mbh objects. 

Given the discussion above, it seems rather peculiar that a 
simple straight line provides a good fit to the Rev 40 variability 
vs. Mbh relation, even when we use 80 and 40 ks long segments 
and for objects with Mq less than 10^ Mq. In principle, we would 
expect a higher spread around the best fit line at low Mbh- There 
are 5 objects with Mbh less than 10^ M© in the 20 ksec Rev re- 
lation. Three of these AGN (namely NGC4395, NGC4051 and 
NGC3227) actually have a low m (ranging from about 0.1 to 2 
%). The other two AGN (MRK766 and NGC4593) are accreting 
at about 60 and 10 % Eddington and, in fact, they lie below the 
best fit relation. Although this is in line with the Model B pre- 
dictions, the presence of only 5 AGN in this Mbh range is ham- 
pering us from deriving any strong conclusion. The CAIXAvar 
(20 ks) has more objects in this Mbh range. As expected, most of 
the objects do lie below the best fit line, suggesting a flattening 
of the relation at low Mbh- We also observe that most of the ob- 
jects with low cTj^i^j, have high m, while the two object with higher 
variability (1H0707-495 and IRAS 13224-3809) have quite un- 
certain Mbh and LboI- Obviously, we need to observe more low 
mass objects (Mbh < 10^), with a high accretion rate, to de- 
termine whether the variability vs. Mbh relation stays linear (in 
log-log space) at low Mbh or not. If it does, our assumptions 
on the PSD need to be modified. On the other hand, if our PSD 
"results" are correct, then the relation should "break" and one 
could not use the best-fit results we give to estimate Mbh below 
the numbers we give above. 




Log (Mg^) Log (Mg^) Log {M^^) Log (Mg^) 

Fig. 16. Left panel: Expected cr^^ vs. Mbh relations for AGN 
with an accretion rate of 0.01 (red lines) and 1 (blue lines) in 
Model B scenario and cr^^^ computed with 250 s bins and 80 
ks intervals. Middle left: Same as before but for 40 ks intervals. 
Middle right: Same as before but for 20 ks intervals. Dashed 
lines show the expected cr^^ computed with time bins of 50 s. 
Right panel: 10 ks intervals. 

8.4. cr~^^ VS. accretion rate 

It is not easy to determine observationaUy if the variability am- 
plitude depends on accretion rate as well. Although reverbera- 
tion mapping can deliver Mbh estimates with uncertainties of a 
factor three, to measure the accretion rate, one requires not only 
a reliable measurement of Mbh, but of LboI as well. However, 
there are larger uncertainties associated with the LboI estimates 
in most objects, which introduce a large scatter in the variability 
vs. accretion rates plots. A second source of significant scatter 
in such plots is the fact that, for a given light curve length, cr^j, 
will be different for different Mbh objects. Consequently, the 
variability vs. ih relation will also depend on Mbh- 

Early studies akeady suggested a weak dependence of cr^^^ 
on accretion rate (Bian & Zhao 2003; Lu & Yu 2001). O'Neill 
et al. (2005) showed that even this weak dependence disappears 
when cr^^j, is "normalized" to Mbh- Zhou et al. (2010) confirmed 
the O'Neill et al. (2005) result and they found that to reproduce 
the small scatter in the cr^^^ vs. Mbh relation, the slope of the 
variability - accretion rate relation has to be smaller than 0.2. 
We found a strong, almost linear correlation between cr^,^^ and 
accretion rate in CAIXAvar However, this is mainly introduced 
by the large Mbh range of the AGN present in this sample, as 
clearly suggested by the fact that, in agreement with the previ- 
ous studies, this correlation is disappeared when we multiply the 
normalised excess variance values with the Mbh of the objects. 

This result is at odds with the McHardy et al. (2006) and 
Koerding et al. (2007) scaling laws for the PSD break time scale, 
if the PSD amplitude was the same in all objects. We found that 
the same PSD model which can explain the small scatter in the 
variability vs. Mbh relations (Model B), also predicts a much 
shallower variability vs. accretion rate relation. If the PSD am- 
plitude decreases with m as PSDamp oc m^" **, and the PSD break 
frequencies scale with Mbh and m according to McHardy et al. 
(2006) and Koerding et al. (2007), then ct^^j^^xMbh °^ m"'^, in 
agreement with the Zhou et al. (2010) results, but also with our 
data (solid lines in the bottom panel of Fig.|9ll 

Clearly, it is difficult to detect such a shallow relation. As 
we discussed above, the most obvious difficulty is the scatter 
introduced in the data plots due to the large uncertainties of 
the accretion rate measurements when we go beyond the 10-15 
nearby better studied AGN. Another important effect is that, in 
data plots with objects which have significantly different Mbh, 
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even the "normalized" estimates ct^^^xMsh are not expected to 
strictly follow the same relation at all m (an effect that was not 
realised by O'Neill et al. 2005, and Zhou et al 2010). For a given 
light curve length of say T, and for sufficiently small Mbh ob- 
jects, we expect that Vbi- > 1/T above a certain accretion rate. 
At higher rates, the cr^j,xMBH vs. accretion rate relation is not 
linear any more (in the log-log space), as the blue and red, solid 
lines in the bottom panel of Fig. |9] clearly demonstrate. This ef- 
fect, is expected to increase even further the scatter in these plots, 
rendering even more difficult the detection of the intrinsic vari- 
ability - accretion rate relation in AGN. 

In any case, the Model B curves shown in both panels of 
Fig. |9]indicate that our results are consistent with the McHardy 
et al. (2006) and Koerding et al. (2007) scaling laws, but only if 
PSDamp oc m This relation introduces a dependence of vari- 
ability on accretion rate, which although weak, is necessary to 
explain other observational results, such as the variability prop- 
erties of the NLS 1 . Previous work on large samples of NLS 1 do 
find a significant difference between the variability of broad and 
narrow line objects that can not be justified by the difference in 
their Mbh, only (Nikolajuk et al. 2009). We also observe indi- 
cations for higher variability amplitude in NLSl, even when we 
consider NLSl and BLSl of comparable Mbh- But we showed 
that, once the difference in their accretion rate is taken into ac- 
count properly, i.e. when we consider the Model B PSDamp de- 
pendence on accretion rate, then both NLSl and BLSl show 
comparable variability amplitudes. This result indicates that the 
variability mechanism is the same in both classes objects, and 
the differences in their variability amplitudes can be fully under- 
stood once the differences in their average Mbh and accretion 
rate is taken into account properly. 



The AGN in CAIXAvar show a strong correlation between cr^^^ 
and the 2-10 keV spectral index F. Early works on this relations 
(Turner et al. 1999; Bian & Zhao 2003) already found indica- 
tions for these two quantities to be correlated (~ 99 % confi- 
dence), with the more variable NLS 1 having steeper spectral in- 
dexes than broad line AGN. However, larger samples indicated 
this correlation to be rather weak (~ 96 % confidence; O'Neill et 
al. 2005). Using a larger sample than ever before (to study AGN 
X-ray variability on time scales less than a day) we estimate the 
significance of the correlation to be as high as 99.99 % with a 
steep slope of ~ 9.2 + 1.2. 

The physical origin of this correlation is not easy to under- 
stand, but we argued that it is the product of two underlying cor- 
relations: the one between variability and accretion rate, and the 
relation between F and m. The latter correlation has been sug- 
gested by many AGN studies in the past (e.g. Porquet et al. 2004; 
PiconcelU et al. 2005; Shemmer et al. 2006; Saez et al. 2008; 
Papadakis et al. 2009; Sobolewska & Papadakis 2009; Risaliti et 
al. 2009; Scott et al. 2011) and has also been observed in BHB 
(see e.g. Wu & Gu 2008). We note that this correlation is not 
significant in CAIXA (Bianchi et al. 2009) but this may be due 
to the flatness of the correlation and the large scatter of the data 
points (both in accretion rate and F) in CAIXA. If F oc as 
Sobolewska & Papadakis (2009) suggested, then since cr^^^ ~ m 
(see Table 2) in CAIXAvar (a relation which is mainly driven 
by the presence of different Mbh objects in this sample, as we 
discussed in the previous section), we would expect to observe a 
cr^j^j. ~ F'" relation, which is entirely consistent with our obser- 
vations. 



If we multiply cTj^js with Mbh, then the variability vs. spec- 
tral index correlation disappears. But in this case, it is difficult 
to argue that this is due to the fact that both cr^^^ and F corre- 
late with Mbh- In fact, both in CAIXA and in CAIXAvar a F vs. 
Mbh correlation is not significant. Instead, when normalized to 
Mbh, the expected variability vs. accretion rate relation is much 
flatter (as we discussed above) and hence much more difficult 
to measure accurately. Therefore, just like with the cr^^^, xMbh 
vs m plot, the ct^^^xMsh vs F plot is flat, but nevertheless, the 
predicted relations, assuming Model B and the Y - m relation of 
Sobolewska & Papadakis, are entirely consistent with the data 
(see Fig.[T4b. 



9. Conclusions 

We studied the variability properties, computing the excess vari- 
ance, of all the radio quiet. X-ray unobscured, AGN observed 
by XMM-Newton for more than 10 ks in pointed observations. 
The XMM-Newton high sensitivity and long orbit allows us to 
measure the variability amplitude of 161 AGN and probe the 
variability on shorter timescales down to 10 ks. We sample 
timescales going from 80 to 10 ks, with the latter being signif- 
icantly smaller than the time scales considered in all previous 
works (like e.g. O'Neill et al. 2005; Zhou et al. 2010). CAIXAvar 
is the largest sample in which short time-scale variability (on 
time-scales smaller than about a day) has been studied and is 
composed primarily by local AGN. Thus the relations and vari- 
ability properties investigated in this work are representative of 
a population of AGN characteristic of the local universe. 

- We systematically explore the energy dependence of the ex- 
cess variance, comparing the variability of the hard (2-10 
keV) with the soft (0.3-0.7 keV) and medium (0.7-2 keV) 
bands. Surprisingly, although different spectral components, 
such as the soft excess and warm absorbers "dominate" the 
emission at softer energies, the soft band variability ampli- 
tude is very well correlated with the hard band variability. In 
particular, no deviations from the linear one-to-one correla- 
tions are observed, suggesting that these components add a 
minor contribution to the total variability. 

- We observe a highly significant and tight (0.4 dex) corre- 
lation between cr^^^ and Mbh for the reverberation sample. 
A larger scatter (0.7 dex) is observed for CAIXAvar This 
is due to the larger uncertainties on the Mbh of the non- 
reverberation estimates. This result implies that our best-fit 
Rev sample results could be used to measure Mbh from crl^^ 
measurements. Provided that these estimates are as accurate 
as the variability estimates for the objects in our Rev sample, 
the Mbh estimates should be more accurate than the ones 
based on single epoch spectra. Particular care though should 
be given to the results from the use of this method if the Mbh 
turns out to be smaller than ~ 5 x 1 O*" - 5 x 1 0^ M© (depending 
on the length of the light curve in use), as in this regime the 
variability - Mbh relation may have a larger intrinsic scatter. 

- The crl^ vs. Mbh correlation has a small scatter (smaller 
than expected, see i.e. McHardy et al. 2006; Koerding et al. 
2007) because the PSD normalisation decreases with accre- 
tion rate as PSDamp °^ m""^. The combination of this relation 
with the PSD break frequency scaling relations of McHardy 
et al. (2006) and Koerding et al. (2007) results in the PSD 
high frequency integral being weakly dependent on accre- 
tion rate, being thus a good tool to estimate Mbh (Nikolajuk 
et al. 2004; 2006; 2009). 
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We note that the universal constancy of the high frequency 
tail (rather than the break frequency) of the PSD in all BH 
accreting sources, seems to suggest a deep link with some 
common physical process experienced by all accreting BH. 

- We measure Mbh or provide stringent upper limits for 65 
AGN, the majority of which (44) have no either reverberation 
mapping or stellar velocity dispersion Mbh estimates and for 
6 of which we did not find any Mbh estimate in literature. We 
also estimate lower limits to the Mbh for the remaining 96 
AGN. 

- We observe a significant correlation between cr^j and m, 
which dissapears when we "normalize" cr^^ to the Mbh of 
the objects, in agreement with the previous results of O'Neill 
et al. (2005) and Zhou et al. (2010). We argue that, the 
lack of significant intrinsic scatter in the cr^j vs. Mbh rela- 
tion, and the lack of a significant correlation between vari- 
ability (normalised to Mbh) and m are in agreement with 
the hypothesis that: all AGN have the same PSD shape, 
with the break frequencies scaling with Mbh and m as in 
McHardy et al. (2006) and Koerding et al. (2007), but only if 
PSDamp a: m""'^, in a fashion similar to what is observed in 
BH binaries (Axelsson et al. 2005; Gierlinski et al. 2008). 

- We observe that NLSl are more variable than broad line 
AGN. This is due to their smaller Mbh and higher accre- 
tion rate. Once rescaled, no difference from broad line AGN 
is required to explain their different variability properties. 

- We observe significant correlations between cr^^^ and both 
the bolometric luminosity and the FWHM of the H/S. Both 
these correlations disappears when the cr^^^ vs. Mbh is con- 
sidered, suggesting that the formers are just by-products of 
the variability vs. Mbh dependence. 

- We observe, for the first time at very high significance (> 
99.99 %), a correlation between cr^^^ and 2-10 keV spectral 
index. This correlation provides indirect support to previous 
suggestions that the spectral index in AGN depends on accre- 
tion rate. In particular, the CAIXAvar cr^j - F relation can 
be fully explained if F oc np-^, as suggested by Sobolewska 
& Papadakis (2009). 

In summary, all our results are in agreement with a picture 
where, to first approximation, all AGN have the same variabil- 
ity properties once rescaled properly for Mbh and accretion rate. 
This confirms and extends the work of McHardy et al. (2006) 
and Koerding et al. (2007) who studied the PSD of less than two 
dozen sources. In our case, this picture is based on the study of 
a significantly larger number of objects. Our work demonstrates 
the significance of studies which are based on the estimation of 
the excess variance, which although as a statistic carries far less 
information than a full PSD analysis, its use allows the consid- 
eration of large AGN samples, hence establishing the validity of 
the current paradigm for the X-ray variabiUty of AGN, and even 
allowing the unexpected detection of important "complications", 
such as the PSD amplitude dependence on accretion rate. 
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Table 5. List of all best fit relations of the CAIXAvar sample as well as their probabilities. For the correlations with significance higher than 65 % we also report the best fit slope 
and normalisation as well as the associated l-cr errors. We do not report the best fit values of the less significant correlations because of their loose physical meaning. 
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central limit theorem, the distribution of the mean of all the ex- 
cess variance estimates will become normally distributed for a 
large number of measurements. Obviously it is not easy to es- 
timate how "large" the number of cr^^^ estimates should be but, 
following this line of reasoning, in the case of sources with more 
than 10 light curve intervals, we consider the error on the mean 
to be symmetric (as if its distribution were normal), and in order 
to estimate the 90 % confidence level, this value is then multi- 
plied by 1.6. Despite the numerous limitations, at least this ap- 
proach allows us to estimate the uncertainty of our final estimate 
without the need to adopt any assumptions about the underlying 
variability process. 



Fig. A.l. Distribution of the observed excess variances com- 
puted for 10 ks intervals for 3C273. The distribution of each 
of the other 6 sources with longer XMM-Newton exposures are 
displayed on the on-line version of the paper The points above 
each distribution indicate the average excess variance and the 
"error bars" around this mean indicate the mean upper and lower 
90 % confidence limits on cr^^^, estimated as described in the 
text. They are compatible with the region where ~ 90 % of the 
points are located (see text for more details). 

Appendix A: Estimating tlie excess variance 
uncertainty 

The uncertainty on the excess variance depends both on the mea- 
surement uncertainties (e.g. Poisson noise) in the light curve and 
on the stochastic nature of the variability. As shown by Vaughan 
et al. (2003) through MonteCarlo simulations, or Ponti et al. 
(2004) in the particular case of large number of photons, the for- 
mer can be approximated by the formula (see eq. 1 1 of Vaughan 
et al. 2003): 



2 < cr2 >) 



<0-f> 2f v: 



(A.l) 



where < crj > is the mean of the uncertainties squared and Fym 
is the fractional variability. 

The uncertainty owing to the stochastic nature of the pro- 
cess is more difficult to estimate. Vaughan et al. (2003) showed 
that the uncertainty increases with the steepness of the power 
spectrum slope and it is large and highly non-Gaussian for steep 
Power Spectral Densities (PSD). A simple estimate of the uncer- 
tainty cannot be analytically derived, however two approaches 
can be pursued. 

A. 1. Estimating tlie "red noise" scatter: more tlian 10 
intervals 

The first approach is to estimate the uncertainty directly from the 
data, deriving the scatter from the distribution of the measured 
excess variances of the different intervals. Every excess variance 
measurement is an independent variable with identical distribu- 
tion (unless the process is non-stationary). This is only approx- 
imately valid, because successive light curve parts may not be 
truly independent, if the "memory" of the system is longer than 
the segment's length. Nevertheless, every time that more than 10 
valid segments are available for the same source, we estimate the 
uncertainty on the mean excess variance computing the sigma of 
the mean of the different measurements. Furthermore, due to the 



A.2. Estimating tlie "red noise" scatter: 1 interval 

When the number of valid segments is just one, we estimated 
the uncertainty on cr^^^ using the results from the Monte Carlo 
simulations of Vaughan et al (2003) and assuming a PSD shape 
for each source. Detailed studies of PSD in AGN have shown 
an almost ubiquitous PSD shape characterised by a steep power 
law shape (a ~ 2) above a special frequency Vbi and a flatter 
(of ~ 1) power law slope below (Uttley & McHardy 2005; but 
see also the case of ARK564; McHardy et al. 2007). In particular 
it has been observed that the characteristic break frequency (vbr) 
scales primarily with mass and with accretion rate (McHardy et 
al. 2006). We thus expect that the slope of the PSD, within the 
frequency range probed by our hght curves (from 4 x 10"^ to 
either 10 ^ 5 x 10"^ or 2.5 x 10"^ Hz, for the 10, 20 and 40 ks 
segments, respectively), may be different for the sources in the 
sample, depending on their Mbh and accretion rate. Assuming 
that all AGN show the same PSD and that the PSD shape scales 
with mass and accretion rate as measured by McHardy et al. 
(2006), we can predict the position of the break frequency and 
the slope of the PSD in the frequency band which corresponds 
to 250 s binned light curves, of duration equal to 10, 20, 40 and 
80 ks. 

Using then the results listed in Table 1 of Vaughan et al. 
(2003), we assume a Mog{S^) = H-0.45 and -0.71, for the pos- 
itive and negative error on cr^^^ , respectively, if the break time- 
scale is longer than the length of the interval (this means that 
in the frequency window on which the excess variance is com- 
puted, the PSD has a slope of -2). We assumed Alog(S^) - +0.28 
and -0.36 if the break time-scale is shorter than the light curve 
time bin (in which case we would expect a PSD of a -1 slope 
in the frequency range sampled by each light curve segment). If 
the break time-scale falls within the frequency window, then we 
combine the errors with the formula: 0.28 x 



Log(v )-Log(v ) i'^^'^ similarly for the negative error). 
When no Mbh estimates are available, we are conservative and 
assume the largest uncertainties, associated with a PSD slope of 
-2. As computed by Vaughan et al. (2003) this scatter estimates 
the 90 % confidence interval. In this way we are able to conser- 
vatively estimate the scatter in the excess variance measurements 
introduced by the red noise, in the case when we have just a sin- 
gle excess variance value. 



A.3. Estimating the "red noise" scatter: checking the 
approximations 

In order to judge how accurate are our excess variance uncer- 
tainties, given the various assumptions that underlie the methods 
we described above, we performed the following test. We con- 
sidered the 7 sources with the longest XMM-Newton observa- 
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tions in CAIXAvar. The number of 10 ksec segments for each 
one of these objects is significantly larger than 10 (so the 10 
ksec excess variance measurement we list in Table |2] for these 
objects is the mean of all the individual measurements, and its 
error is based on the true scatter of the points around their mean, 
as explained above). Figure lA.ll shows the distribution of the 
observed 10 ksec excess variance values, for these objects. As 
expected, whenever we have a "signal" (i.e. the cr^^^^ values are 
positive) the distributions are asymmetric with a tail at large val- 
ues. In the case of 3C273, due to its large Mbh, we do not expect 
large amplitude variations on such short time scales. As a result, 
the distribution of the lOks cr^^^^ values is dominated by the un- 
certainties associated with the Poisson noise of the light curves, 
hence it is more symmetric. 

We then used the method described above to obtain the un- 
certainty of each individual cr^^^ value, and we calculated the 
mean upper and lower 90% confidence limits on cr^^^^, for each 
object. The points on top of each sample distribution in Fig. lA.ll 
indicate the average excess variance, and the "error bar" around 
this mean indicate the mean upper and lower 90% confidence 
limits on cr^^^. A visual inspection shows that these confidence 
limits are compatible with the area where ~ 90% of the points 
are. In fact, we measure that the average 90 % uncertainty es- 
timated with our method actually contains 84, 92, 80, 93, 87, 
97, 92 % of the measured values of NGC4051, 1H0707-495, 
MGC-6-30-15, MRK766, MRK335, NGC3516 and 3C273, re- 
spectively. Thus the typical difference between the estimated 90 
% uncertainty and the one measured from the observed distribu- 
tion is less than 10 %. 

The 1-0" "error" has been estimated simply dividing the 90 
% uncertainties by 1.6. This approximation is valid for normal 
distributions, only. However, in the cases of NGC4051, 1H0707- 
495, MGC-6-30-15, MRK766, MRK335, NGC3516 and 3C273, 
using the observed distribution, we can compare the 1-cr uncer- 
tainty computed in this way with the corresponding probabilities 
that are 73, 69, 51, 62, 58, 54 and 82 %, respectively. We, thus, 
measure that the typical difference between the 1 -cr computed in 
this way and the 1-cr measured from the observed distribution, is 
of the order of 10-15 %, reaching maximum values of 25 % for 
MCG-6-30-15. 

Fig. lA.ll suggests that the method we have adopted to es- 
timate the uncertainty on the excess variance in the case when 
there is just a single interval available results in an acceptable- 
conservative estimate of the true scatter of cr^^^. The same figure 
also suggests that this result should be valid for AGN with differ- 
ent Mbh and intrinsic excess variance values which spans almost 
four orders of magnitude. 



Appendix B: o-^^^ computed in fixed iengtli 
intervals and time dilation with redshift 

Since CAIXAvar is composed mainly of local AGN (see Fig.[Tll, 
the impact of the differences in red-shift on the excess variance, 
if we estimate it using intervals with a fixed length in the ob- 
servers frame, should be minimal. This is even more obvious in 
the case of the sources with at least one variable segment de- 
tected, whose redshift is less than 0.3 in all (but one) cases. For 
these objects, segments of fixed length in the observer's frame 
should imply differences of less than ~ 30% in the rest frame 
segment's length. Since the excess variance does depend on the 
maximum frequency sampled, we expect that the intrinsic excess 
variance should be different, depending on the source's redshift. 
However we verified that, for power-spectra like the ones typi- 
cally observed in AGN (see discussion in Section 6.1.3 below), 
the resulting differences in the intrinsic variance of the sources 
should be less than ~ 20% (for a large range in BH mass and ac- 
cretion rates). This maximum difference is much smaller than the 
observed scatter in all the variability plots we study below. For 
this reason, we decided to work with the excess variance mea- 
surements that we estimated from the intervals we mentioned 
above, whose length is fixed, irrespective of the sources redshift. 
Regarding the highest redshift source (at z=0.9) in the sample 
of "variable" objects, we always checked that its presence does 
not affect in anyway our results from the study of the correla- 
tions between the variability amplitude and the various physical 
parameters we present in this paper. 



A.4. Estimating tlie "red noise" scatter: between 2 and 10 
intervals 

When the number of valid segments is higher than one and lower 
than 10, we estimated the stochastic scatter for each segment as 
detailed above, in the case of sources with just a single interval 
available. Then, following O'Neill et al. (2005), our final esti- 
mate of the mean excess variance is equal to the square root of 
the sum of the squared "error" of the individual segments, di- 
vided by the number of intervals. Finally, both in the case of 
single and less than 10 intervals, we combined in quadrature 
the stochastic scatter and the one associated with the Poissonian 
noise (see eq. lA.U . 
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